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1*.  SUmaNCNTARV  NOTATION 


IX  W2  uJT+M( 


- %»  Remits  ara  raportad  fron  tha  first  yaar  of  a  three-year  Joint  BRL/AMMRC/SRI  program 

to  develop  a  computational  capability  to  pradlct  tha  behind -t ha-armor  fragment  anvironmant 
for  spaced  armor  attacked  by  long  rod  penetrators.  Tha  materials  chosen  for  study  ware 
depleted  uranium-  and  tungst da-alloy  penetrators,  rolled  homogeneous  armor,  and  an 
alec troslag -remalt -treated  ateel  armor.  Phenomenology  experiments  indicated  that  adiabatic 
shear  banding  is  the  dominant  microstructural  failure  mods  underlying  target  plugging  and 
fragmentation  aa  well  ae  penetrator  nose  erosion.  A  previously  developed  computational 
model  for  shear  banding  way  improved,  calibrated  with  dynamic  material  property 
experimental  data,  and  applied  to  preliminary  computational  simulations  of  normal  impact 
experiments.  In  addition,  an  approximate  formula  was  developed  to  estimate  the  critical 
strain  for  onset  of  shear  banding.  Promising  agreement  between  these  preliminary 
computations  and  experimental  observations  was  obtained,  Information  was  also  obtained 
on  tensile  failure  of  rolled  homogeneous  armor  by  ductile  void  activity,  and  a  method  - - - 
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rot  generating  high  strain  cat*  plastic  7 laid  lafoxaatloa  was  developed.  Future  work 
will  (attend  thin  approach  to  obllqua  Inpacta,  aultlple  plates,  sad.  If  the  results  of 
ongoing  phenoa apology  aatparinanta  establish  their  Importance,  to  other  alcrostructurel 
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I  INTRODUCTION  AND  SUMMARY 


The  objective  of  this  three-year  joint  3RL/AMMRC/SRI  program  is 
to  develop  improved  material  models  for  computationally  simulating  long 
rod  penetration  of  spaced  armor,  thereby  providing  the  capability  to 
predict  the  behind-the-armor  fragment  environment  for  given  penetrator 
attacks  as  well  as  to  aid  in  armor  and  penetrator  design. 

Our  approach  to  meet  this  objective  has  six  main  parts.  First, 
materials  were  selected  for  study,  and  penetration  phenomenology  experi¬ 
ments  are  being  performed  in  which  rods  are  fired  at  various  angles  of 
incidence  against  target  plates.  The  rod  and  target  specimens  are 
recovered,  sectioned,  and  metallographically  examined  to  reveal  the  key 
microstructural  damage  mechanisms  that  govern  failure  of  both  penetrator 
and  target. 

The  baseline  materials  chosen  for  study  are  depleted  uranium  (DU)- 
3/4  wtZ  titanium  alloy  rods  fired  against  rolled  homogeneous  armor  (RHA) 
targets.  Additional  materials  to  be  studied  are  tungsten  alloy  rods 
and  electroslag-remelt  (ESR)  steel  armor. 

In  the  second  part  of  the  program,  laboratory  material  property 
characterization  experiments  are  being  performed,  where  necessary,  to 
measure  the  properties  that  govern  the  damage  mechanisms.  '  ^ 

Third,  previously  developed  computational  models  of  the  various 
damage  processes  are  being  modified  and  adapted  to  describe  the  observed 
failure  mechanisms. 

Fourth,  these  models  are  calibrated  against  the  experimental  data 
with  computaticnc'  simulations  of  the  material  characterization 
experiments. 


la  the  fifth  part  of  the  approach,  Cha  calibrated  aaterial  models 
are  being  incorporated  in  a  two  dimensional  wave  propagation  code, 

HOP,  developed  by  Dr.  Mark  Wilkine  and  coworkera  at  Lawrence  Livermore 
Laboratory  (LLL).  This  code  aust  be  mod if led  somewhat  to  accomodate 
the  material  models. 

Finally,  in  the  sixth  part  of  the  program  the  ability  of  the  models 
to  correctly  predict  the  main  features  of  long  rod  penetration  of  spaced 
armor  is  being  assessed  and,  if  the  results  are  promising,  the  models 
will  then  be  used  to  predict  behind- the-armor  fragment  environments  and 
to  aid  in  armor  and  penetrator  design. 

The  progress  made  on  the  program  during  the  first  year  is  summarised 
below.  RHA  material  was  obtained  from  3RL  early  in  the  program, 
and  the  metallurgical  characterization  has  been  completed  (see  Tables 
1  and  2,  and  Figures  1-6) .  Quasi-static  tensile  tests  and  dynamic 
expanding  cylinder  testa  have  also  been  performed  by  SRI  on  this 
material,  and  plate  slap  experiments  have  been  performed  by  BRL  personnel. 
Phenomenology  experiments  Involving  single,  monolithic  RHA  targets 
impacted  normally  by  4360  steel  rods  have  been  compel tad  and  analyzed 
(see  Figures  7-9),  and  in  addition,  metallurgical  examinations  have 
been  completed  on  SO  fragments  from  rods  that  have  impacted  spaced  RHA 
targets  at  oblique  angles  of  incidence  (Figures  10  and  11) . 

The  phenomenology  experiments  performed  to  date  indicate  that  the 
primary  micro structural  damage  mode  for  erosion  of  the  penetrator  nose 
and  for  plugging  of  the  target  is  adiabatic  shear  banding.  The  primary 
mlcrostructural  damage  mode  for  rear  surface  armor  damage  (spall  and/or 
petal-ling)  is  probably  ductile  void  nucleatlon,  growth,  and  coalescence. 
The  primary  damage  mode  for  failure  of  the  rear  portion  of  the  DU  rod 
is  as  yet  undetermined,  although  flash  x-ray  pictures  suggest  brittle 
fracture. 

Significant  progress  was  made  this  year  on  adapting  a  previously 
developed  ductile  void  growth  model  to  the  armor  penetration  problem. 
However,  because  of  the  observed  importance  of  adiabatic  shear  banding 
in  the  failure  of  both  penetrators  and  armor,  we  have  concentrated  on 
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•hear  band  nodal  ■edification  and  adaptation  during  the  first  year  and 
are  witholding  judgeaent  on  the  laportance  of  ductile  void  kinetics 
until  all  the  phenomenology  results  are  in. 

Froa  available  data  on  expanding  cylinders  of  4340  steel  (Figures 
12*16) ,  we  confirmed  that  the  onset  of  adiabatic  shear  banding  can  be 
treated  as  a  • lapis  plastic  instability,  in  agreeaent  with  suggestions 
aade  by  aany  previous  workers.  The  instability  is  triggered  when  the 
decrease  lu  yield  strength  caused  by  adiabatic  thermal  softening  in  the 
prelocalisation  notarial  exceeds  the  Increase  caused  by  the  work 
hardening  (Figures  B.l  and  B.2).  The  material  becones  globally  con- 
stltutively  unstable  at  the  equivalent  plastic  strain  at  which  th'  true 
effective  stress  begins  to  decrease  with  Increasing  strain.  Thereafter, 
any  perturbation  will  localise  and  grow.  We  derived  a  simple  formula 
for  roughly  estini  .ng  this  critical  threshold  strain  based  on  the 
yield  stress,  work  hardening  exponent,  and  the  thermal  softening 
coefficient.  This  formula  was  found  to  correctly  rank  the  propensity 
of  various  materials  to  fora  adiabatic  shear  bands  (see  Table  B.l),  and 
we  expect  the  formula  to  be  useful  for  back-of-the-envelope  assessment 
of  armor  and  penetrator  materials. 

Once  the  threshold  conditions  for  shear  band  nucleatlon  are 
fulfilled,  the  computational  model  must  describe  the  shear  band  numbers, 
sizes,  and  degree  of  coalescence  to  predict  the  fragmentation.  The 
associated  nucleatlon,  growth,  and  coalescence  models  have  been  cali¬ 
brated  by  analysis  of  available  data  from  contained  fragmenting  cylincer 
(CFC)  tests  on  4340  steels  of  several  hardnesses  (see  Figures  20-24) . 

The  resulting  models  will  be  applied  during  the  next  year  to  similar 
data  now  being  generated  for  RHA  steel  and  to  be  generated  for  DU. 

Finally,  preliminary  computational  simulations  for  normal  impacts 
of  long  rods  on  single  target  plates  have  been  performed  with  a  simple 
HEMP-llke  wave  propagation  code  (called  TROTT) .  In  these  calculations 
the  improved  adiabatic  shear  band  model  (SHEAR3)  was  tested,  and  the 
feasibility  of  using  SHEAR3  to  activate  preexisting  slide  lines  was 
examined.  These  runs  were  successful  (see  Figures  29-31) .  Furthermore, 
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penetration  mechanisms  such  as  penetrator  bulging  and  impact  of  the 
bulge  shoulder  with  the  target  surface  were  studied.  The  results  showed 
that  double  slide  lines  will  be  required  in  HEMP  to  handle  the  case  of 
extreme  bulging  of  whe  penetrator  r.ose  in  normal  impact,  as  is  expected 
at  impacts  near  the  ballistic  limit,  and  to  handle  the  case  of  oblique 
Impact  when  the  side  of  the  rod  impacts  the  target  surface.  Fortunately, 
Mark  Wilkins  and  coworkers  at  LLL  have  recently  inserted  a  double  slide 
line  algorithm  in  one  version  of  HEMP.  We  have  requested  this  version 
of  HEMP,  and  we  will  implement  the  double  slide  line  algorithm  during 
the  next  year. 

In  the  remainder  of  this  report  we  discuss  the  above  topics  in 
some  detail,  although  frequent  reference  will  be  made  to  References 
1  through  11,  the  relevant  monthly  progress  reports.  In  some  cases, 
for  convenience  of  access,  key  portions  of  these  references  will  be 
-repeated  almost  verbatim.  In  Section  II  we  describe  both  the  phenomenology 
and  material  characterization  experiments.  Section  III  covers  the 
computational  model  development,  the  calibration  of  the  models  with 
experimental  data  (in  particular,  the  calibration  of  SHEAR!  with  4340 
steel  data),  and  the  preliminary  computational  simulations  of  penetration 
experiments.  In  Section  IV  we  discuss  our  conclusions  and  outline  plans 
for  the  coming  year. 


II  EXPERIMENTS 


A.  Material  Selection.  Procureme nt.  and  Metallurgical  Analysis 

Tha  baseline  asCsrisls  seise ted  for  study  in  cooperation  with  BRL 
sod  AMMRC  personnel  were  s  rolled  homogeneous  armor  steel  (RHa)  and  a 
depleted  uraniue  alloy  (DU).  In  addition,  tungsten  alloy  penetrator 
materiel  and  ESR  steel  armor  aatarial  will  be  studied. 

Early  in  the  program,  three  18  inch  x  19  inch  x  4  inch  thick  plates 
of  RHA  material  were  shipped  to  SRI  by  BRL.  The  mechanical  property 
and  metallurgical  characterizations  of  this  RHA  material  are  described 
in  detail  in  Raferances  S  and  8.  The  composition  was  assumed  to  be 
the  same  as  that  measured  from  a  4-inch-thick  plate  of  the  same  heat 
of  RHA  by  Bench22  and  Hauver13,  and  is  given  in  Table  1. 

To  our  knowledge,  the  thermal  properties  of  the  RHA  have  not  been 

measured.  However,  its  thermal  properties  should  be  similar  to  those 
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of  steels  alloyed  in  the  same  manner.  The  Metals  Reference  Book 
reports  a  thermal  conductivity  of  0.082  cal/hr  cm  *C  and  a  specific 
heat  of  0.118  cal/g'C  for  a  3Z  Nl-Cr  low  alloy  steel. 

The  alpha  phase  Hugoniot  of  this  heat  of  RHA  has  been  measured  by 
Hauvar23.  In  plotting  the  wave  velocity  (v)  versus  the  particle 
velocity  (u),  the  Hugoniot  is  given  by:  v  ■  4.51  +  1.43  u.  The 
Hugoniot  of  RHA  is  very  similar  to  the  Hugoniots  of  pure  iron  and  AISI 

steel. 

Metallogrsphlc  analysis3  shoved  that  the  RHA  from  the  4-inch  plate 
is  quite  clean  for  a  commercial  heat  of  such  dimensions.  It  should  be 
rated  at  no  more  than  Field  No.  1  on  the  ASTM  inclusion  content  scale. 

In  general,  the  physical  and  micros true tural  characteristics  of  this 
RHA  appear  reasonably  uniform  and  isotropic  within  a  7.5-cm  (3-inch) 
central  core. 
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Table  1 


COMPOSITION  OF  AN  RHA  PLATE  10  CM  (4  IN.)  THICK* 


Klaaent 

VtZ 

Carbon 

0.27 

Manganese 

0.27 

Phosphorous 

0.001 

Sulfur 

0.008 

Silicon 

3.15 

Nickel 

3.47 

Copper 

0.05 

Chroad.ua 

1.37 

i 

I 

0.01 

Molybdenum 

0.10/0.25 

Aluminum 

0.03 

Iron 

Balance 

i  i2  13 

Taken  from  reporta  by  Benck  and  Hauver. 
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Nevertheless,  the  mechanical  properties  exhibit  significant 

anisotropy  that  is  related  to  ttrt  inclusion  orientations.  To  investigate 

this  anisotropy,  as  well  as  to  study  aicrostructural  fracture  mechanisms 

and  to  confirm  that  our  material  la  essentially  identical  to  that  studied 
12  13 

by  Bench  and  Hauver  ,  we  performed  several  smooth  round-bar  tensile 
tests  on  specimens  machined  in  different  orientations.  These  experiments 
and  the  results  are  discussed  in  detail  in  kafcrence  8;  in  the  following 
paragraphs  we  summer ised  the  procedures  and  the  hey  results. 

The  round-bar  ’’ensile  specimens  were  0.250  inch  (6.35  am)  in 
diameter,  3.750  inches  (9.52  cm)  in  length,  and  with  a  1.0-inch  (25.4 
am)  gage  length.  They  were  machined  from  the  RHA  plate  in  the  three 
principal  orientations:  in  the  plana  of  rolling  and  parallel  to  the 
rolling  direction  ("rolling"  direction) ,  in  the  plane  of  rolling  and 
perpendicular  to  the  rolling  direction  ("long  transverse"  direction), 
and  the  through-the-thichness  direction  ("short  transverse"  direction) . 
The  rolling  direction  was  identified  during  the  previous  metallurgical 
analysis.  The  specimens  were  machined  from  the  plate  starting  1  inch 
(25.4  mm)  in  from  any  surface  to  eliminate  edge  effects  such  as 
hardening  from  the  cutting  torch  or  surface  hardening  from  the  rolling 
process. 

Quasi-static  uniaxial  tensile  tests  were  performed  on  the  specimens 

in  accordance  with  ASTM  Specification  E-8.  In  all  cases  the  specimens 

were  pulled  to  failure  at  a  crosshead  speed  of  0.02  inch/min  (0.51  mm/ 

min)  using  a  standard  Instron  testing  machine.  This  gave  a  nominal 
-4  -1 

strain  rate  of  10  sec  .  All  tests  were  performed  at  room  temperature. 
For  one  specimen  in  each  orientation,  the  plastic  strain  to  failure 
was  measured  with  a  calibrated,  1-inch  (25.4-mm)  gage  length  extensometer 
On  the  remaining  specimens,  axial  and  lateral  plastic  strain  in  the 
neck  region  was  measured  using  single-frame  photography.  From  these 
data,  we  obtained  load-deflection  curves  and  true  stress-true  strain 
curves  for  each  orientation. 
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One  fractured  spec  lam  in  each  orlantaclon  was  sectioned  along  Its 
axia  to  atudp  tha  void  growth  characteristics  in  tha  nack  region  and 
beyond .  In  addition,  the  fracture  surface  of  one  speclaen  in  each 
orientation  waa  exealned  using  scanning  electron  aicroscopy  (SEM) . 

Table  2  suaaarlses  the  test  result*  for  each  specimen.  The  true 
fracture  stress  (TPS)  was  determined  by  dividing  the  load  at  failure 
by  the  cross-sectional  area  of  the  neck  at  failure.  It  is 

apparent  that  the  rolling  and  long  traverse  directions  exhibit  similar 
properties  but  that  the  short  transverse  direction  exhibits  significantly 
degraded  properties,  particularly  in  the  ductility-dependent  properties 
(true  fracture  stress,  elongation  to  failure,  and  reduction  in  area) . 

This  type  of  anisotropy  between  the  rolling  plane  and  tha  through-the- 
thlckness  plane  la  characteristic  of  thick,  aa-rolled  plates. 

Figure  1  shows  the  true  stress-true  strain  behavior  of  specimens 
K-l  (rolling  direction),  L-2  (long  transverse  direction),  and  S-2 
(short  transverse  direction)  as  measured  from  the  photographs  of  the 
neck.  The  loading  path  in  true  stress-true  strain  space  is  same  in  all 
three  directions.  The  only  significant  discrepancy  between  the  rolling 
and  long  transverse  directions  and  short  transverse  direction  is  in 
true  strain  to  failure.  Hence,  the  deformation  in  all  three  directions 
should  be  descrlbable  by  the  same  flow  equation.  These  data  (at  least 
up  to  a  true  strain  of  0.8)  are  descrlbable  by  a  power  lav  hardening 
equation  of  the  fora: 


where  0  -  true  stress,  K  ■  200  ksi  (1.38  GPa),  e  -  true  strain,  and 
N  -  0.19. 

Also  plotted  in  Figure  1  are  the  true  stress-true  strain  quasi- 

statlc  tensile  data  on  4-inch-thick  (10.16-cm)  RHA  plate  reported  by 
12  * 

Benck  at  BKL.  Although  there  is  some  scatter,  the  SRI  data  and  the 

* 

The  values  of  K  and  N  in  Equation  (II-l)  were  calculated  usi,.",  both 
the  SRI  data  and  the  SRL  data. 
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SUMMIT  Or  QUASI -STATIC  TENSILE  TESTS  ON  ENA 


tnnAmmuttim  n>iwvMWMuinnuu  IVJIMrwLMnAMU'MUlf 


E:  Young's  Modulus  TFS :  true  fracture  stress 

Y:  0.22  offset  yield  strength  e-  elongation  to  failure 

UTS:  ultimate  tensile  strength  R.A.:  reduction  in  area 


QUASI-STATIC  TENSILE  TESTS 


itU  collected  at  ML  in  tha  chraa  plats  oriaatationa  agree  raaaonably 
wall  (eoapata  also  Tabla  II  in  Benck*2  with  Table  2  In  thla  raport) . 

The  data  of  ML  also  tad  lea tt  aaiaotroplc  behavior  in  tha  through-the- 
thickness  dlractloo.  It  thua  appears  probabla  that  tha  DBA  plata  tested 
by  SKI  la  from  tha  same  4-inch-thick  (10.16-cm)  plata  of  KHA  taatad  by 
laock  aad  that  spot-to-epot  variations  la  thla  plata  do  not  appaar  to 
ba  excessive.  Therefore,  tha  results  of  dynamic  and  static  failure  mods 
experiments  conducted  by  SKI  on  BHA  can  ba  fairly  compared  with  tha 
reusita  of  tha  experiments  conducted  at  BKL. 

These  tensile  tast  results  indicate  that  this  4-inch-thick  (10.16 
cm)  plata  of  HA  exhibits  significantly  reduced  ductility  in  the  through- 
t ha- thickness  direction  relative  to  the  rolling  plane  directions. 
Moreover,  nlcrostructural  characterisation  indicated  that  this  BRA 
plate  uas  reasonably  isotropic  with  respect  to  grain  also,  grain  shape, 
and  constituent  phases;  bunco,  these  nlcrostructural  characteristics 
cannot  ba  the  cause  of  the  ductility  anisotropy.  Bouev  or,  the  micro- 
structure  was  observed  to  contain  large,  nonuniformly  distributed  ItaS 
inclusions.  The  inclusions  were  either  roughly  spherical  or  oblong, 

10  to  20  un  in  diameter,  and  concentrated  in  dusters  that  were  elongated 
in  tha  direction  of  rolling. 

Figure  2  is  a  scanning  electron  micrograph  of  tha  fracture  surface 
of  specimen  K-l  (rolling  direction).  The  fracture  surface  exhibits  a 
cup-end-cone  type  of  fracture,  that  is,  tensile  void  nudeatlon,  growth, 
and  shearing  of  the  remaining  ligaments  between  the  voids  in  the  center 
of  the  specimen,  followed  by  shearing  around  the  periphery.  Specimen 
L-2  (long  transverse  direction)  exhibited  a  similar  fracture  surface. 
Figure  3  details  the  void  structure  in  the  center  of  specimens  R-l  and 
L-2.  The  voids  in  both  specimens  era  in  general  smell  ('v-  1  ym) , 
spherical,  and  quite  deep.  The  particles  that  nucleated  these  voids 
could  not  be  identified,  but  their  small  sise  suggests  that  they  were 
nucleated  around  the  fine  (*\»  1  ym) ,  uniformly  dispersed  iron  carbide 
(Fe^C)  particles  characteristic  of  a  balnltic  microstructure. It 
thus  appears  that  the  fracture  process  in  the  rolling  and  long  trans¬ 
verse  sp.cimens  la  controlled  by  the  nucleation  and  growth  of  voids 
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SHOWING  THE  CONVENTIONAL  CUP-AND-CONE  TYPE  OF  F 


(a)  SPECIMEN  R-1  (Rolling  Diraction) 


(b)  SPECIMEN  L-2  (Long  Transvarsa  Diraction) 


MP-7893-28 

FIGURE  3  DETAILS  OF  THE  VOID  STRUCTURE  IN  THE  CENTER 

OF  THE  FRACTURE  SURFACES  OF  SPECIMENS  R-1  AND  L-2 

13 


1 

ibibsocBQininiiDii^^ 


around  tho  Mali,  uniformly  dlaparsod  Fe^C  partlcloa.  Tha  largar  MnS 
lndualona  that  ara  eoncantracad  along  planaa  parallal  to  tha  rolling 
dir act Ion  do  not  appaar  to  algnlflcantly  affact  tha  fracture  process  of 
tha  rolling  and  long  transveraa  apaclaana. 

Plgura  4  la  a  scanning  electron  micrograph  of  tha  fracture  surface 
of  specimen  8-2  (abort  transverse  direction) .  Approximately  60X  of  tha 
fracture  surface  exhibits  a  ductile  dimple  characteristic,  and  approxi¬ 
mately  40X  exhibits  a  shear  lip.  Tha  aaynmetric  appearance  of  the 
fracture  surface  suggests  that  the  ductile  dimple  fracture  process  was 
controlled  by  the  nucleation  and  growth  of  voids  around  a  nonuniform 
distribution  of  large  inclusions.  Figure  5(a)  shows  some  of  the  large, 
elongated  dimples  found  on  the  flat  part  of  the  fracture  surface.  The 
particles  visible  at  the  bottom  of  the  voids  were  Identified  as  MnS  by 
nondlspersive  x-ray  analysis.  The  large  dimples  visible  in  Figure  8(a) 
are  larger  than  any  found  on  the  fracture  surfaces  of  the  rolling  plane 
specimens  (Figure  3).  However,  the  fracture  surface  of  the  sheer  lip, 
Figure  5(b),  Is  very  similar  in  appearance  to  the  fracture  surfaces  of 
the  shear  lips  on  the  rolling  plane  specimens  (R-l  and  L-2) .  It  thus 
appears  that  the  lower  ductility  observed  in  tho  short  transverse 
specimens  can  be  attributed  to  the  nucleation  and  growth  of  large  voids 
across  at  least  part  of  the  neck  region,  followed  by  plastic  tearing  of 
any  remaining  ligament.  It  also  appears  that  these  large  voids  were 
nucleated  around  large  MnS  particles  observed  to  be  concentrated  in 
dusters  elongated  in  the  direction  of  rolling. 

Figure  6  shows  the  polished  cross  sections  of  one  fracture  specimen 
in  each  of  the  three  principal  directions.  The  sire  and  distribution 
of  voids  in  the  vicinity  of  the  neck  in  specimens  R-l  (Figure  6(a)] 
and  L-2  [Figure  6(b)]  are  quite  similar.  As  expected,  the  void  size 
and  distribution  are  much  different  in  specimen  S-2  [Figure  6(c)]  in 
which  a  number  of  large  voids  are  concentrated  primarily  under  the  flat 
portion  of  tha  fracture  surface.  It  is  these  voids  that  reduce  the 
ductility  of  the  short  transverse  specimens  relative  to  the  rolling 
and  long  transverse  specimens.  Note  also  that  the  voids  in  specimens 
R-l  and  L-2  are  elongated  in  a  direction  parallel  to  the  load,  while 
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FIGURE  4  SEM  OF  THE  FRACTURu  SURFACE  OF  SPECIMEN  S-2  (Short  Transverse  Direction) 

Approximately  60%  of  the  fracture  surface  exhibits  a  ductile  dimple  characteristic,  and 
approximately  40%  is  covered  by  a  shear  lip  (compare  with  Figure  2). 


DETAILS  OF  THE  FRACTURE  SURFACE  OF  SPECIMEN  S-2  (Short  Transverse  Direction) 

(a)  Shows  the  large  elongated  ductile  dimples  characteristic  of  60%  of  the  fracture  surface. 

(b)  Shows  that  the  shear  lip  is  covered  by  small  shallow  dimples  distorted  by  shear. 


the  voids  in  ipcclam  S-2  srs  slongatsd  in  a  diraccion  perpendicular 
to  tha  load.  This  indicates  that  tha  large  MaS  inclusions  in  tha  short 
Cransvarsa  spac loans  ara  elongated  in  a  direction  perpendicular  to  tha 
load  axis  (parallel  to  Che  rolling  direction)  and  that  tha  voids  nucleated 
by  such  inclusions  had  not  grown  ouch  by  plastic  flow  (had  not  elongated 
in  the.  loading  direction)  at  tha  tine  of  fracture. 

Va  conclude  that  the  observed  aechanical  anisotropy  is  due  to 
clusters  of  large  MnS  inclusions  elongated  in  tha  rolling  direction. 

The  fracture  process  in  the  rolling  and  long  transverse  directions  is 
controlled  by  the  nucleation  and  growth  of  voids  around  tha  fine 
(<  1  ya)  unlfornly  dispersed  Fe^C  particles  characteristic  of  bainltic 
micros true tur as,  followed  by  shearing  of  the  remaining  ligaments.  The 
fracture  process  in  the  short  transverse  direction  is  controlled  by  the 
nucleation  and  growth  of  voids  around  the  large  MnS  inclusions  elongated 
in  the  direction  of  rolling,  followed  by  shearing  of  the  remaining 
ligaments. 

B.  Phenomenology  Experiments 

1.  Response  of  RHA  to  Normal  Impacts  of  Long  Steel  Bods 

We  performed  a  series  of  experiments  in  which  flat-nosed  4340 
rods  (Rc50)  were  fired  at  normal  incidence  into  RHA  plates.  The  rods 
were  1/4  inch  (0.635  cm)  thick  and  1/2  inch  (1.27  cm)  long.  The  RHA 
target  plates  were  machined  from  the  plate  discussed  above  and  were 
3  inches  (1.18  cm)  square  and  1/4  inch  (0.635  cm)  thick.  The  Impact 
velocities  were  chosen  to  range  from  just  below  to  just  above  the 
ballistic  limit,  and  the  targets  were  recovered  and  sectioned  after  the 
shots  for  metallographic  examination. 

The  results  are  shown  in  Figures  7  through  9.  Several 
observations  can  be  made  immediately.  First,  the  RHA  material  is 
fairly  resistant  to  shear  bandirg.  Nevertheless,  the  target  failure 
mechanism  is  shear  banding  and  resultant  plugging.  Second,  shear 
banding  is  the  erosion  mechanism  for  the  nose  of  the  4340  steel 
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Projectile  Note 
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MP-7M3-32A 

FIGURE  7  POLISHED  AND  ETCHED  CROSS  SECTION  THROUGH  CENTER  OF 
RHA  TARGET  IMPACTED  BY  A  BLUNT-NOSED  4340  STEEL  PRO¬ 
JECTILE  AT  2150  FT/SEC  (655  m/sec) 

Not*  abaanc*  of  shear  banding  in  targat  material. 


(b) 


1/4" 


Shear 

Banded 

Penetrator 

Nose 

Fragments 
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FIGURE  8  IMPACT  SIDE  OF  (a)  AND  POLISHED  AND  ETCHED  CROSS  SECTION 
(b)  THROUGH  CENTER  OF  RHA  TARGET  IMPACTED  BY  BLUNT- 
NOSED  4340  STEEL  PROJECTILE  AT  2330  FT/SEC  (710  m/sec) 
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1/4” 


Embedded 

Penetrator 


1/4" 


(b) 


MP-7893-42A 


FIGURE  9  BLUNT-NOSED  4340  STEEL  PROJECTILE  WITHIN  PLUG  ON  REAR 
SURFACE  OF  TARGET  PLATE  (a)  AND  POLISHED  AND  ETCHED 
CROSS  SECTION  THROUGH  RHA  TARGET  AND  PLUG  FOLLOWING 
IMPACT  AT  2600  FT/SEC  (792  m/sec)  (b) 


penetrator a.  Node*  that  the  mushrooming  nose  fragments  causa  a 
dlacontlnuous  or  scalloped  crater  In  the  target.  Finally,  the  mush¬ 
rooming  effect  of  the  fragmenting  penetrator  noae  combined  with  the 
realatance  of  the  IHA  to  ahear  banding  and  plugging  has  caused  the 
onset  of  shear  banding  in  the  RHA  to  occur  first  at  a  radius  around 
the  impact  center  that  is  significantly  greater  than  the  original  radius 
of  the  penetrator.  (Figure  9  ahovs  a  case  where  the  plug  radius  is 
about  33Z  larger  than  the  original  penetrator  radlua.) 

2.  Response  of  RHA  to  Oblique  Impacts  of  Long  DU  Rods 

In  cooperation  with  BRL  and  AMflLC  personnel,  we  have  procured 
DO  penetrator  roda  for  use  in  penetration  phenomenology  experiments 
involving  Impacts  at  varioua  velocities  and  angles  of  Incidence  against 
single  and  multiple  BHA  target  platea.  These  DU  roda  are  described  in 
Section  V.  Posttest  recovery  and  metallographic  examination  of  sectioned 
target  and  penetrator  fragments  will  be  made  to  study  the  microstructural 
failure  modes.  A  key  consideration  ia  the  Influence  of  the  failure 
rechanlans  on  scaling.  For  example,  brittle  fracture  of  the  rear  portion 
of  the  penetrator  using  bending  stresses  is  expected  to  be  sensitive  to 
the  plastic  sone  slse, 

«c  "  A(Klc/0y)2  *  (II~2> 

wh-'re  A  is  a  geometry-dependent  coefficient,  1L  is  the  fracture 
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toughness,  and  oy  is  the  yield  stress. 

At  the  date  of  this  writing,  these  phenomenology  experiments 
have  not  yet  been  completed.  However,  preliminary  information  has  been  ) 

obtained  from  metallographic  examination  of  samples  from  experiments  « 

performed  at  BRL  in  which  DU  penetrators  were  fired  obliquely  at  RHA 
targets.  Figure  10  shows  polished  sections  of  the  rear  portion  of  a 
65-g  DU  rod  that  had  penetrated  scaled  triple  armor  with  an  impact 
velocity  of  985  m/sec  and  a  residual  velocity  of  562  m/sec.  (The  angle 
of  Incidence  is  unknown) .  Figure  11  shows  a  section  through  a  RHA 

i 
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target  in  which  a  DO  fragment  is  embedded.  (Tbs  inpact  velocity  and 
tha  angle  of  incidence  are  unknown.) 

It  la  dear  froa  Figures  10  and  11  that  adiabatic  shear 
banding  was  the  key  failure  aechanisa  of  the  penetrator  nose  in  both 
canes.  Of  particular  interest  is  the  band  shown  in  Figure  11.  Here 
the  aartensltlc  platelets  are  seen  to  bend  sharply  as  they  enter  the 
bend,  ehowing  clearly  the  localisation  of  plastic  atraln. 

The  shear  bands  in  Figures  10  and  11  were  also  seen  to  contain 
vo Ida,  la  egresment  with  observations  reported  by  Roger s.*^ 

In  summary,  although  definite  conclusions  nust  await  completion 
of  the  phenomenology  experiments,  it  appears  that  shear  banding  is  the 
key  damage  mechanism  for  both  the  RHA  target  and  the  nose  of  the  DU 
penetrator.  However,  we  expect  that  ductile  void  growth  will  also  be 
Important  for  high  impact  velocity  cases  for  which  back  surface  spall 
occur. 


3*  Response  of  RHA  to  Impacts  of  Long  Tungsten  Alloy  Rods 

AMMRC  personnel  are  planning  to  perform  a  series  of  oblique 
impact  experiments  with  long  tungsten  alloy  rods.  Discussions  have 
been  held  with  John  Meacall  of  AMfRC  to  determine  if  posttest  recovery 
of  fragments  will  be  possible  to  examine  failure  phenomenology.  In 
addition,  BRL  personnel  plan  to  perform  normal  impact  experiments  with 
tungsten  alloy  rods  (L/D  ■  10)  fired  against  2-inch-thick  RHA  target 
plates  that  are  Instrumented  with  stress  gages.  Discussions  with 
George  Hauver  of  BRL  have  been  held  to  plan  for  posttest  recovery  of 
penetrator  and  target  fragments  for  failure  phenomenology  examinations 
at  SRI. 

Thus,  tungsten  alloy  penetrator  failure  phenomenology  data 
should  be  forthcoming  during  the  next  year. 
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(b»  View  of  the  microetructuro  of  the  ponotrator 
shaft  back  from  the  impact  and.  Tha  hotoa 
moat  likafy  raauitad  from  oxide  inchmons 
strung  out  along  tho  extrusion  axis  and  than 
pulled  out  during  grinding  and  oloctropolishing. 

MP-7893-70 

FIGURE  10  RECOVERED  URANIUM-3/4  wt%  TITANIUM  AFTER  PERFORATION 

OF  SCALED  TRIPLE  ARMOR  (First  Plata:  0.24  cm;  Second  Plate:  0.68  cm 
Third  Plate:  1.26  cm) 

The  impact  velocity  was  985  m/sec  and  the  residual  velocity  was  562  m/$ec. 
The  residual  mass  was  12.8  g. 


FIGURE  11  STEEL  TARGET  PLATE  No.  4  IMPACTED  AT  AN  OBLIQUE  ANGLE  BY  A  URANIUM-3/4  wt% 
TITANIUM  ALLOY  PENETRATOR 


C.  Material  Characterisation  tggrtjgtg 


1.  Contained  Fragment ins  Cylinder  (CTC)  Experlme nti 

Th«  CTC  wptrlunti  war*  developed  In  *n  **rll*r  project*® 

Co  auiuri  th*  ptf«Mt«ra  In  our  computational  nodal  for  ahaar  banding. 
Tha  technique,  daacrlbad  in  datall  in  Kafaranca  18,  an tall*  datonatlng 
explosive  in  a  thick-walled  tuba  of  tha  aatarial  of  Intaraat.  Th* 
tuba  ia  aurroundad  by  concentric  layer*  of  plaatlc  and  ataal  respectively, 
so  that  th*  expansion  of  tha  inner  tuba  can  ba  halted  at  various  atagas 
of  development  by  varying  tha  thickness  of  tha  plaatlc  layer.  The 
explosive  density  can  also  be  varied,  and  th*  strain  histories  can 
thereby  be  tailored  by  varying  tha  explosive  density  and  th*  plastic 
layer  thickness.  Tha  strain  histories  throughout  tha  specimen  tub* 
are  calculated  with  our  two-dimensional  wave  prepagatlon  code,  TROTT. 

To  date,  four  CFG  experiments  have  been  performed  using  RHA 
specimen  tubes  having  a  length  of  3-3/4  Inches,  an'l.D.  of  1-1/8  inches, 
and  a  wall  thickness  of  1/4  inch.  Th*  experimental  parameters  and  some 
qualitative  results  are  shown  in  Table  3.  Th*  recovered  specimen  tubes 
wars  sliced  in  half  along  a  plan*  containing  the  axis,  and  the  final 
wall  thickness  was  measured  as  a  function  of  axial  position.  The 
equivalent  plaatlc  strain  (ep)  was  calculated  by  th*  formula  eP  -  /3* 
e?  (1  +gr)*  where  cr  (th*  radial  plastic  strain)  equals  the  natural 
logarlthmrof  the  ratio  of  th*  final  wall  thickness  to  th*  original  wall 
thickness  and  c  (axial  plastic  strain)  equals  the  fractional  increase 
in  length  of  th*  specimen.  Figure  12  shows  the  results  for  all  four 
specimen  tubes  in  terms  of  equivalent  plastic  strain  versus  axial 
position. 

Metallogrephic  examinations  have  been  completed  on  th*  RHA 
specimen  tubes  L2  and  L4 ,  whose  axes  were  parallel  to  the  long  trans¬ 
verse  direction  of  the  83A  slab.  Figure  13  contains  two  views  of  one 
half  of  the  specimen  L4,  showing  large  shear  band  damage  extending 
almost  th*  entire  length  of  th*  tube,  but  concentrated  in  two  specific 
circumferential  regions,  and  lesser  damage  adjacent  to  these  regions. 

A  polished  end  etched  section  perpendicular  to  the  axis  of  the  tube 
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MA-7893-45A 


FIGURE  12  EQUIVALENT  PLASTIC  STRAIN  DETERMINED  FROM  MEASUREMENTS 
OF  RHA  SPECIMEN  TUBES  RECOVERED  FROM  CFC  EXPERIMENTS 
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FIGURE  13 


TWO  VIEWS  OF  HAi_F  OF  THE  RHA  SPECIMEN  TUBE  L4  RECOVERED  FROM 
CFC  EXPERIMENT 

(Plane  in  Figure  14  is  cut  at  a  distance  of  2-3/8"  from  bottom  of  specimen  tube.) 
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at  a  dlatanca  of  2-3/8  Inches  from  Che  botCoa  end  is  shown  in  Figure  15. 
Except  for  a  few  saall  bands,  all  of  Che  shear  bands — froa  the  incipient 
bands  barely  noticeable  along  the  inside  radius  to  the  couple  of  large 
bands  that  have  extended  completely  through  the  tube — have  a  direction 
that  Is  approximately  parallel  (within  about  15*)  to  the  rolling 
direction.  This  anisotropy  is  caused  by  planes  of  shear  weakness  in 
the  rolling  direction  due  to  the  rolling  process. 

We  then  machined  the  specimen  tubes  labeled  T1  and  T3  with 
the  tube  axla  in  the  short  transverse  direction,  as  shown  in  Figure  15. 
It  is  clear  from  Figure  15  that  shear  bands  that  cut  across  the  rolling 
direction  should  be  easier  to  produce  in  the  short  transverse  specimens, 
because  in  these  specimens  the  plane  of  weakness  Induced  by  rolling  will 
be  closer  to  the  mode  2  orientation,  as  defined  in  Figure  16,  and  the 
tube  expansion  in  the  mode  2  orientation  produces  only  about  half  the 
resolved  shear  strain  produced  in  the  mode  1  orientation. 

Indeed,  the  recovered  specimens  labeled  T1  and  T3  in  Table  3 
showed  a  large  number  of  mode  1  incipient  bands  distributed  uniformly 
around  the  circumference;  the  bands  appeared  only  after  equivalent 
plastic  strains  greater  than  about  80Z,  as  compared  to  a  nudeation 
threshold  strain  of  less  than  50Z  for  the  bands  in  the  rolling  plane 
in  the  L  specimens. 

The  following  conclusions  can  therefore  be  drawn  from  work 
completed  to  date:  (1)  there  is  a  significant  anisotropy  in  the 
resistance  of  RHA  to  shear  banding;  (2)  bands  in  planes  parallel  to 
the  rolling  plane  nucleate  at  equivalent  plastic  strain  (ep)  less  than 
50Z,  whereas  bands  in  planes  perpendicular  to  the  rolling  plane  nucleate 
(£P)  80Z;  and  (3)  our  computational  shear-banding  model  may  have  to 

account  for  this  anisotropy.  t 

In  armor  penetration  (specifically  by  a  blunt-nosed  projectile 
impacting  at  right  angles)  plugging  is  caused  by  growth  of  shear  bands 
in  planes  perpendicular  to  the  rolling  plane.  These  bands  can  be 
activated  in  CFC  specimens  with  axes  parallel  to  the  short  transverse 
direction  of  the  RHA  slab  (T-specimen) .  Incipient  shear-band  levels 


30 


FIGURE  15  CFC  SPECIMENS  MACHINED  FROM  RHA  SLAB,  SHOWING  MODE  1  SHEAR 
BANDS  IN  PLANES  PERPENDICULAR  TO  ROLLING  PLANE 


ORIENTATION  MODE  1 
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Orientation  mode  2 
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ORIENTATION  MODE  3 
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FIGURE  16  GEOMETRY  AND  NOMENCLATURE  FOR  SHEAR  BANOS  IN  CONTAINED- 
FRAGMENTING-CYLINOER  EXPERIMENTS 


have  been  achieved  thus  far;  shots  are  nov  under  way  to  produce  higher 

daaage  levels. 

Shear  banding  in  planes  parallel  to  the  rolling  plane  may  be 
important  in  other  armor-penetrator  geometries  (e.g.,  oblique  impact 
of  hemispherical -nosed  or  ogival -nosed  projectiles),  particularly  in 
the  case  of  RHA,  which  is  less  resistant  to  shear  banding  in  those 
planes.  These  bends  can  be  activated  in  CFC  specimens  tilth  axes 
parallel  to  the  long  transverse  direction  (L-apeclmen) .  If  future 
phenomenology  experiments  show  that  these  bands  are  important,  additional 
CPC  experiments  will  be  performed  with  L-speclmens  to  obtain  incipient 
to  moderate  levels  of  damage  (higher  levels  of  damage  have  already 
been  obtained)'. 

2.  Sv—etric  Taylor  Testa 

As  will  be  discussed  in  Section  III,  the  dynamic  (adiabatic) 
stress-strain  curve  above  yield  la  an  important  component  in  our  com¬ 
putational  damage  models.  The  onset  of  shear  banding  is  particularly 
sensitive  to  the  dynamic  work  hardening  and  adiabatic  thermal  softening. 
In  addition,  ductile  void  growth  and  the  size  of  the  plastic  zone 
(which  governs  the  resistance  of  the  penetrator  rod  and  target  plate 
to  brittle  fracture;  see  Eq.  11-2)  are  strong  functions  of  the  dynamic 
yield  function.  We  therefore  need  to  obtain  yield  data  at  the  rates  of 
Interest  for  RHA  and  DU  (as  well  as  for  tungsten  alloy  and  ESR  steel 
later  in  the  program) . 

The  strain  rates  of  Interest  range  from  10^  sec  *  to  10^  sec 

a  relatively  inaccessible  aud  largely  unexplored  region  of  measurement. 

The  two  experimental  methods  that  appear  to  hold  the  most  promise  for 

19 

exploring  this  strain  rata  region  are  the  dynamic  torsion  test  and  a 

20  21 

modification  of  the  Taylor  test.  ’ 

The  dynamic  torsion  test  is  currently  being  developed  and 

adapted  to  measurements  on  armor  and  penetrator  materials  under  con- 

22 

tract  with  BRL  by  Minneapolis  Honeywell  Corporation.  Therefore,  in 

20 

our  work  we  have  concentrated  on  a  modified  version  of  the  Taylor  test. 
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The  classical  Taylor  tsst  involves  impacting  a  rod  on  a  rigid  surface. 

Our  modified  version  entails  the  symmetric  impact  of  two  identical 
right  circular  cylinders,  each  about  3/8  inch  in  diameter  and  1*1/2 
inch  long  (thaae  dimensions  may  be  varied) .  Figure  17  shows  a  schematic 
drawing  of  the  experisMnt.  One  of  the  cylinders  is  accelerated  by  a 
light  gas  gun  to  a  velocity  in  the  range  of  0.5  to  1.0  km/aec,  where 
it  impacts  a  stationary  cylinder.  Both  cylinders  are  supported  by 
stiff,  but  weak  fixtures,  so  that  the  supporting  fixtures  interface 
negligibly  with  the  wave  propagation  and  deformation  resulting  from 
the  impact.  Immediately  following  impact,  a  high-speed  framing  camera 
records  a  series  of  photographs  at  approximately  1,000,000  frames/sec 
to  produce  a  record  of  the  "mushrooming"  plastic  deformation  at  the 
Interface  of  the  two  projectilea.  Posttest  recovery  examination  of 
projectiles.  Including  metallographic  examination  of  any  internal  shear 
banding  or  tensile  failure,  completes  the  data.  The  experiment  is  then 
computationally  simulated  in  Iterative  calculations  with  a  two-dimensional 
computer  code.  The  dynamic  yield  surface  parameters  are  varied  until 
the  experimental  and  computational  "mushrooming"  agree. 

The  symmetric  Taylor  test  has  two  advantages  over  the  classical 
Taylor  test.  The  first  is  the  elimination  of  possible  computer  simulation 
problems  regarding  the  movement  of  the  cylinder-target  interface  and 
any  possible  yielding  or  failure  in  the  actually  nonrlgid  target.  The 
second  advantage  is  the  elimination  of  the  uncertainty  of  whether  the 
projectile  sticks  to  or  slides  freely  along  the  target  interface,  a 
problem  that  again  leads  to  computational  uncertainities.  In  the 
syvmetrlc  Taylor  test,  the  interface  between  the  two  cylinders  does 
not  move  axially  in  ths  center-of-mass  coordinate  system,  and  the  two 
cylinders  stick  precisely  together  as  their  interface  expands  rapidally 
outward.  We  note  that  to  attain  the  same  impact  conditions  in  the 
synnetrlc  Taylor  test  must  be  approximately  twice  as  high  as  that  in 
the  classical  version.  However,  this  is  not  a  problem  for  the  metals 
of  Interest  in  this  program. 
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Tc  data,  we  have  par fora ad  two  symmetric  Taylor  tasts  to  test 
the  experimental  procedure.  Figure  18  shows  the  framing  camera  records 
and  the  recovered  epaclaana  from  a  test  with  aluminum  rods.  The  experl- 
mental  parameters  ware  chosen  to  closely  simulate  a  classical  Taylor 
test  performed  by  Wilkins  at  al.  on  the  same  material.  After  the 
teat  the  rods  ware  recovered,  sectioned,  and  examined  to  sea  if  the 
mushrooming  was  due  completely  to  homogeneous  plastic  flow,  as  desired, 
or  whether  heterogeneous  failure  mechanisms  had  been  active.  Figure  19 
shows  that,  for  this  experiment,  the  petaanent  deformation  was  due 
entirely  to  homogeneous  plastic  flow,  as  expected  from  Wilkins's  work. 

In  future  work  this  technique  will  be  applied  to  RHA  and  DU 
rods.  In  preparation  for  the  analysis,  Brian  Scott  of  BRL  plans  to 
test  our  procedure  for  obtaining  the  dynamic  yield  function  by  running 
HEMP  computer  simulations  of  the  above  aluminum  rod  Impact  experiment 
and  comparing  the  results  with  the  framing  camera  data. 

3.  Smooth  Bound  Bar  Tension  Tests 

The  smooth  round  bar  tension  tests  on  XHA  material  were 
discussed  earlier  in  Section  IIA  in  connection  with  the  mlcrostructural 
tensile  failure  mechenism,  namely,  anisotropic  nudeatlon  of  voids  at 
MnS  inclusions,  followed  by  growth  and  coalescence  to  cause  fracture. 
These  tests  also  give  us  the  quasi-static  yield  function  to  use  as  a 
starting  point  for  the  Iterative  calculations  of  the  symmetric  Taylor 
tests  for  evaluation  of  the  dynamic  yield  function. 

Similar  round  bar  tension  tests  of  DU  material  say  be  performed 
in  the  future  if  available  data  are  found  to  be  ambiguous. 

♦.  Plate  Impact  Tests 

Plate  impact  tests  designed  to  produce  short-lived  tensile 

pulses  in  the  target  material  under  conditions  of  uniaxial  strain  are 

useful  for  studying  dynamic  tensile  failure  phenomenology.' 

Dr.  Gerald  Moss  has  performed  a  series  of  such  experiments  at  BRL  with 
24  25 

tungsten  alloy  and  RHA  material,  and  he  is  designing  similar 
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FIGURE  19  PHOTOGRAPHS  OF  RECOVERED  ALUMINUM  ROD  SHOWING  THAT 
THE  PERMANENT  DEFORMATION  WAS  DUE  ENTIRELY  TO  HOMO¬ 
GENEOUS  PLASTIC  FLOW 


npttiMti  with  DU  to  bt  performed  if  the  phenomenology  experiments 
rnul  ductile  void  growth  to  bo  on  iaportent  failure  mechanism  for  the 
DO  ponotretore. 

The  rooulto  fro*  thooo  experiments  will  bo  analysed  to  provide 
input  to  the  MAC  nodal  for  dynamic  ductile  void  growth. 


Ill  ANALYSIS 


A.  Computational  Modal  Adaptation 

The  computational  material  modal  used  in  tha  armor  penetration 
simulations  auat  have  the  potential  to  deacribe  elastic-plastic 
deformation,  quasi  static  and  dynamic  tensile  fracture,  and  adiabatic 
shear  banding.  In  the  following  paragraphs  we  describe  these  portions 
of  the  model. 

1.  Elastic-Plastic  Model 

The  elastic-plastic  modal  we  are  now  using  applies  a  von  Mises 
yield  function  with  linear  work  hardening.  Nonlinear  work-hardening 
rate-dependent  models  are  available  if  needed.  No  new  development  of 
elastic-plastic  models  was  undertaken  during  this  program;  instead, 
we  have  applied  these  models  currently  available  in  our  wave  propagation 
computer  codes. 

2.  Ductile  Fracture  Model  for  Intermediate  Strain  Rates 

During  projectile  penetration,  fracture  can  occur  at  several 
positions  along  the  penetrator  axis  and  on  the  back  surface  of  the 
target.  In  all  cases  the  fracture  probably  occurs  in  10  to  100  usee. 
For  usual  penetrator  and  target  materials,  the  fracture  occurs  by  a 
ductile  mode  and  is  initiated  by  the  growth  of  voids,  although  brittle 
fracture  may  also  occur.  We  have  begun  to  develop  a  ductile  fracture 
model  that  will  permit  detailed  simulation  of  these  damage  processes. 

The  model  includes: 

*  Nudeation  or  activation  of  new  voids  either  as  a  function 
of  stress  and  time  (as  in  DFRACT26)  or  as  a  function  of 
plastic  strain  (as  in  DFRACTS22) . 

*  Elastic  void  expansion  and  contraction,  as  outlined  by 
Seaman  et  al.2§ 


Void  compaction  under  compression,  es  described  by 
Seemen  at  al.28 

A  static  fracture  threshold  as  described  by  Gurson.2^ 


With  these  features,  the  model  can  simulate  static  and 
dynamic  ductile  fracture  over  the  complete  range  of  strain  rates  con¬ 
sidered  in  this  study. 


The  major  effort  was  devoted  to  the  static  fracture  threshold, 
which  was  the  suet  difficult  factor  to  incorporate.  The  threshold 
problem  is  discussed  below,  and  an  outline  is  presented  of  a  technique 
for  combining  the  four  items  listed  above  into  a  single  model. 


To  resolve  the  threshold  problem,  we  used  a  simplified  form 
of  Gurson's  model2^.  Gurson  treated  a  sphere  of  rigid-plastic  material 
containing  a  spherical  void.  As  the  tenalle  stress  state  was  increased, 
regions  of  plastic  deformation  formed  around  the  void.  The  void  growth 
threshold  was  reached  when  the  plastic  tones  had  extended  to  the  point 
that  void  expansion  could  occur  without  resistance  by  the  remaining 
rigid  tutorial.  Stress  states  from  pure  shear  to  uniform  tension 
(hydrostatic)  were  examined.  The  threshold  stress  function  was  written 
as  a  yield  curve  where 
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Gurson  also  shows  that  a  normality  flow  rule  describes  the 
stress-strain  relations  for  the  material  in  the  yielding  range.  The 
normality  rule  is  a  condition  between  the  stresses  (average  values 
over  the  boundary  of  the  sphere  of  material)  on  the  yield  curve  and 

p 

the  average  plastic  strains  throughout  the  sphere. 

rd0ijd£1j  "  0  (HI-2) 


This  normality  flow  rule  is  commonly  used  in  some  form  in  plasticity 
problems . 

In  our  use  of  Gurson* s  model,  we  simplified  the  yield  function 
to  the  ellipse  as  follows: 


o2  9  P2 

Yz (1  -  f)2  +  4Y^ (in  f)2  "  1 


(III-3) 


This  ellipse  has  the  correct  intercepts  and  normality  properties  at 
the  Intercepts. 


The  requirements  of  the  normality  rule  for  an  isotropic 

material  governed  only  by  the  invariants  O  and  P  lead  to  a  relation 

between  the  plastic  strains.  The  ratio  S  of  the  stress  invariant 
P  — P 

increments  dP  and  da  corresponding  to  these  plastic  strain  increments 
is 
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(III-4) 


where  $  is  the  yield  function,  G  and  K  are  the  shear  and  bulk  moduli, 
— N 

respectively,  a  and  P  are  invariants  computed  from  the  total  strain 
Increments,  and  H  and  P  are  computed  from  the  elastic  components  of 
the  strain  Increments.  The  usual  normal  to  the  yield  function  $  has 
the  slope 
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By  comparing  Eqa .  (III-4)  and  (III-5),  w*  saa  chat  tha  slopa  corre¬ 
sponding  Co  Cha  normality  rule  Is  proportional  to  tha  normal  slopa  of 
tha  yield  curve.  With  tha  normality  rule  for  tha  plastic  stress-strain 
relations,  a  solution  procedure  for  stresses  under  incremental  strain 
loading  can  be  formulated. 

Tha  solution  procedure  for  stress  caused  by  any  strain  incre¬ 
ment  la  developed  by  first  assembling  the  conditions  that  must  be  met. 
These  are  reduced  to  a  single  implicit  equation  for  pressure.  The  so¬ 
lution  is  then  obtained  by  a  regula  falsi  technique.  The  four  conditions 
that  ssist  be  satisfied  simultaneously  are: 

(1)  Tha  plastic  strains  in  the  strain  interval  must  be 
normal  to  the  mean  point  on  the  yield  surface  during 
tha  interval. 

(2)  Tha  stress  point  must  be  on  the  yield  surface. 

(3)  Tha  void  volume  is  determined  by  the  pressure  and 
internal  energy  in  the  solid  material. 

(4)  Tha  yield  stress  is  determined  by  tha  cumulative  plastic 
strain. 


These  four  requirements  determine  O,  P ,  f,  and  Y.  Below,  each  require¬ 
ment  is  derived  in  analytical  form  and  then  combined  to  form  a  single 
equation  for  pressure  P. 


The  normality  condition  outlined  above  leads  to  a  relation 
between  the  plastic  components  of  O  and  P  and  the  stress  normals  Sq  and 
Sj  at  the  initial  and  final  states  of  the  interval. 
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where  Ac  and  AP  are  plastic  components  of  o  and  P  .  Expressions  for 
SQ  and  are  obtained  from  Eq.  (III-4)  and  the  function  $  in  Eq.  (III-3) 
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E  is  internal  energy.  According  to  Carroll  and  Holt^®,  the  average 
pressure  P  on  the  cross  section  is  reduced  front  Pg  by  the  ratio  of 
gross  to  solid  density 


(III-ll) 


Solving  Eqs.  (III-10)  and  (III-ll)  for  Vg  ■  l/pg,  we  obtain 

1  -  PVp  /K 

Vg  -  - (III-12) 

po  -  VE/KS 

The  void  fraction  is  then 


»  -  PC0/KS 


f0  -  <VE/Ks 


(in-13) 


This  expression  for  void  fraction  naturally  accounts  for  elastic  and 
thermal  expansions  of  the  void  because  both  effects  are  present  in  the 
solid  volume  calculation. 


The  expression  for  pressure  P  can  now  be  written  by  substi¬ 
tuting  cr  from  Eq.  (III-8)  into  Eq.  (III-3).  The  result  is 


where  Y  and  f  are  obtained  from  the  associated  Eqs.  (III-9)  and  (III-13) , 
and  the  value  of  P  in  the  radical  is  from  a  previous  estimate. 

The  solution  procedure  then  has  the  following  steps: 

(1)  Estimate  an  initial  pressure  P  . 

A 

(2)  Compute  the  void  fraction  f  from  Eq.  (III-13). 

(3)  Compute  A,  B,  and  C. 
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(4)  Compute  pressure  P  from  Eq.  (III-14). 

(5)  Reestimate  Pft  and  return  to  step  2  if  P  and  P  are  not 
sufficiently  close . 

If  ongoing  phenomenology  experiments  Indicate  quasi-static 
▼old  growth  to  ba  an  Important  damage  mode,  we  will  then  Incorporate 
this  ductile  fracture  model  into  a  computer  subroutine  during  the 
next  year. 

3.  Fracture  Models  for  High  Strain  Rates 

23 

In  previous  work,  dynamic  fracture  models  were  developed 
that  describe  the  nudeation,  growth,  and  coalescence  of  statistical 
numbers  of  brittle  cleavage  cracks  or  ductile  plastically  growing  . 
voids.  In  this  program  we  have  not  developed  these  models  further, 
but  will  adapt  them  in  their  current  form  by  calibrating  them  with 
the  data  obtained  from  the  BRL  plate  impact  requirements. 

4.  Adiabatic  Shear  Banding  Model 

As  discussed  above,  the  phenomenology  experiments  performed 
to  date  indicate  that  adiabatic  shear  banding  is  the  single  most  impor¬ 
tant  microstructural  damage  mechanism  for  long  rod  penetration  because 
it  governs  both  the  erosion  of  the  penetrator  nose  and  the  plugging  of 

the  target.  Therefore,  we  have  devoted  a  large  effort  to  further 

10 

developing  and  calibrating  our  previous  shear  band  model.  The  pre¬ 
vious  shear  band  model  subroutine,  discussed  in  Reference  18,  is  named 
SHEAR2,  whereas  the  new,  upgraded  subroutine  is  named  SHEAR3. 

The  following  changes  were  made  to  convert  SHEAR2  to  SHEAR3 : 

■  The  shear  band  size  distributions  were  changed  from  a 
series  of  discrete  points  to  an  analytical  form  (see 
Reference  27  and  Appendix  A). 

■  The  stress  relaxation  associated  with  damage  was  com¬ 
pletely  rederived  (see  Reference  4  and  Appendix  A). 

*  The  nucleatlon  and  growth  laws  were  generalized  to  allow 
more  strain  rate  dependence  (see  Reference  11  and  Section 
III-B) . 


•  The  number  of  special  variables  per  cell  required  by  che 
subroutine  was  reduced  (see  Reference  11  and  Appendix  A). 

•  The  threshold  conditions  for  initiation  of  shear  band 
activity  were  derived  from  plastic  instability  theory 
and  a  simple  foraula  was  derived  for  predicting  the 
critical  strain  for  shear  banding. 

The  above  changes  are  described  in  a owe  detail  in  Appendices 
A  and  B.  Appendix  A  describes  the  SHEAR3  subroutine  and  Appendix  B 
describee  the  derivation  of  the  threshold  conditions  for  onset  of  shear 
banding. 


A  key  result  discussed  in  Appendix  B  is  the  critical  strain 
for  onset  of  adiabatic  shaar  banding,  as  given  by  a  foraula  of  the  form 


(111-15) 


where  ?  is  the  critical  equivalent  plastic  strain;  and  n  is  the 
hardening  exponent  in  a  hardening  law  of  the  form 


T  -  T  (e*)®  ,  (III-16) 

o 

where  T  is  the  effective  shear  yielding  stress;  En  is  the  specific  melt 
energy;  p  is  the  density;  and  a  is  a  linear  thermal  softening 
coefficient  froa  the  relation 


T  "  1  "  T  »  (III-17) 

a 

where  V  is  the  specific  plastic  strain  energy  density. 

The  Eq.  (III-15)  clearly  shows  the  competitive  effects  of 
strain  hardening  and  theraal  softening,  and  expresses  the  criterion 
that  the  shear  banding  will  begin  as  a  plastic  instability  when  the 
global  thermal  softening  caused  by  adiabatic  plastic  deformation 
overwhelms  the  global  work  hardening.  Appendix  B  also  gives  a  formula 
for  more  complicated  thermal  softening  behavior.  Finally,  Appendix  B 
compares  predicted  critical  strains  for  onset  of  shear  banding  for 
several  metals  with  observations,  and  shows  that  the  correct  ranking  is 
predicted.  * 
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B.  Generation  and  Tastlni  of  Shear  Banding  Modal  Parameters 


The  material  parameters  for  the  shear  banding  nodal ,  SHEAR3, 
ara  ganaratad  at  praaant  almost  anciraly  froa  containad  f ragman ting 
cyllndar  (CFC)  data.  In  tha  currant  staga  of  tha  program,  va  hava 
data  baaaa  for  two  aatarlala,  4340  ataal  and  RHA  ataal.  CFC  data 
for  DO  will  ba  ganaratad  In  futura  work.  By  far  tha  most  axtanslva 
data  baaa  la  that  for  4340  ataal  of  thraa  hardnaasaa;  thasa  data  ware 
ganaratad  by  SRI  in  a  pravloua  program  for  AMMRC^.  In  tha  following 
paragrapha  wa  daacrlba  tha  currant  atata  of  analysis  of  t»ie  data  for 
both  4340  and  RHA  ataal. 

1.  Analyala  of  4340  Data 

Extensive  ahaar  banding  data  wara  available  from  six 
contalned-fragmenting-cy Under  axparimanta  on  4340  ataal  hardened  to 
Bc40.  Wa  hava  taken  aeveral  atapa  to  uaa  4340  data  to  aid  In  the 
development  and  verification  of  tha  ahaar  band  modal.  In  Reference  7 
wa  made  a  first  attempt  to  deduce  initiation,  nudeatlon,  and  growth 
descriptions  for  tha  ahaar  bands.  In  that  effort  wa  assumed  that  the 
ahaar  band  nudeatlon  and  growth  ratea  should  correlate  with  strain 
rata  histories  that  peaked  about  20  uaac.  after  detonation.  However, 
later  examination  of  our  computational  simulations  showed  chat  Che 
deduced  threshold  strains  of  about  0.2  were  not  reached  until  later, 
when  the  strain  rate  was  significantly  lower.  Therefore,  we  next 
reduced  the  shear  band  data  to  a  fora  where  nudeatlon  and  growth 
processes  could  ba  more  dearly  Identified  and  quantified.  The  shear 
band  modal  was  than  modified  to  account  for  these  observed  processes. 
The  detonation  process  In  TROTT^,  the  two-dimensional  wave  propagation 
coda  used  for  tha  fragmenting-round  simulations,  has  also  been  modified 
to  batter  represent  tha  explosion  in  low-density  PETN.  To  avoid 
tangling  of  the  computational  calls  in  the  numerical  representation  of 
tha  explosive,  an  automatic  rezoner  has  been  written  and  debugged. 

The  shear  band  data  reduction,  shear  band  model  changes,  and  the 
detonation  algorithm  are  described  here.  The  rezone  method  will  be 
described  later  after  we  have  gained  more  experience  in  using  it. 
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The  4340  data  vara  obtalnad  by  counting  and  measuring  the 
ahaar  bands  observed  on  the  surfacaa  of  six  cylinders  after  contained- 
fragmenting- round  experiments.  The  experimental  configuration  is 
shown  in  Figure  20.  The  buffer  thickness  and  the  density  of  the  explo¬ 
sive  were  varied  to  provide .a  range  of  strain  rates,  peak  strains,  and 
levels  of  shear  banding. 

Following  counting  of  the  shear  bands,  the  cumulative  numbers 
were  plotted  in  graphs  like  Figure  21.  Then  each  distribution  was 
fitted  to  the  exponential  relation 


Nf  -  Nq  exp  (-R/Rj) 


(111-18) 


where 


is  the  cumulative  number  of  bands /cm 
with  a  radius  greater  than  R 

2 

is  the  total  number  of  bands/cm 

la  a  sixa  parameter  for  the  distribution. 


Vlth  the  distributions  reduced  to  two  characterizing  parameters,  Nq  and 
R^,  plots  can  be  made  to  study  the  Initiation,  nuclaatlon  rate,  and 
growth  rate  of  the  bands. 


The  initiation  and  nucleatlon  rate  of  shear  bands  was  examined 
by  plotting  Nq  as  a  function  of  plastic  strain  for  <11  the  cylinders t 
as  shown  in  Figure  22.  The  plastic  strain  was  obtained  by  measuring  the 
thickness  change  in  the  cylinders  or  large  fragments  after  the  test. 

As  derived  in  Reference  33,  equivalent  strain  is 


iP- 


V  ep 

rijEij 


•  #<*;> 


Px2 


Px2 


♦  <*e> 


+  <e;r) 


( III— 19) 


where  eQ,  are  strains  in  the  radial,  circumferential,  and  axial 
directions,  respectively.  For  the  calculations,  we  assume  that  the 
axial  and  volumetric  plastic  strains  are  zero,  so  that 
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FIGURE  20  CROSS  SECTION  OF  CONTAINED-FRAGMENTING-CYLINDER  EXPERIMENTS 

Buffer  thickness  X  varies  from  0.96  to  1.9  cm,  PETN  density  p  varies  from  1.20  cm  to  1.47  •fan* . 


SURFACE  DENSITY  OF  SANDS  WITH  LENGTH  >  L 


0  0.2  0.4  0.6  0.8  1.0 

SHEAR  BAND  LENGTH.  L  —  cm 

MA-6401-6 

FIGURE  21  SHEAR  BAND  DISTRIBUTION  FOR  SHOT  2 
IN  4340  STEEL.  R{40 


52 


TOTAL  NUMBER  OF  SHEAR  BANOS/cm2  IN  THE  SIZE  DISTRIBUTIONS 


0.22  0.24  0.26  0.28  0.30  0.32  0.34  0.36 

EQUIVALENT  PLASTIC  STRAIN.  ?p 


MA-7M3-41 

FIGURE  22  MEASURED  NUMBERS  OF  SHEAR  BANDS/cm3  VERSUS  PLASTIC  STRAIN 
BASED  ON  WALL  THICKNESS  MEASUREMENTS  FOR  SIX  CONTAINED- 
FRAGMENTING-CYLINDER  EXPERIMENTS  IN  4340  STEEL.  R£40 
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Then 


*p  •  VT  i«;i 


For  a  change  In  thickness  frosi  aQ  to  a^t  the  average  radial  strain  is 


and 


(111-20) 


cP 


(III— 2 1) 


By  Its  definition  fron  the  thickness  change,  c  represents  the  plastic 
strain  taken  homogeneously  by  the  material  and  not  the  strain  that  has 
been  taken  by  the  shear  benda. 


-P 

An  alternative  derivation  of  c  can  be  made  based  on  the 
circumferential  strain 


o 


(III-22) 


where  rQ  and  r^  are  initial  and  final  radii  at  some  point  in  the 
cylinder.  Let  r  be  the  average  radius  of  the  cylinder.  Then  the 
average  circumferential  strain  is 


e 


e 


(III-23) 
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(III-24) 


CP  -  t  VF‘n 

Mciun  this  strain  is  basad  on  tha  notion  of  tha  axtarnal  and  internal 
surfaces,  it  should  represent  the  total  strain  imposed  by  the  explosive 
and  not  merely  the  strain  taken  homogeneously.  Values  of  this  second 

•p 

definition  of  C  were  calculated  in  a  few  cases  and  compared  with  the 
result  obtained  from  Eq.  (III-21).  The  radius-based  strain  always 
exceeded  the  thickness-based  value  by  5Z  or  10Z.  For  the  cylinders 
with  such  heavy  damage  that  fragments  are  formed,  the  radius-based 
strain  cannot  be  obtained  reliably,  so  the  imposed  strain  is  not  known. 
Therefore,  the  thickness-based  strain  was  used  in  the  plots  as  a  rough 
measure  of  the  imposed  strain. 

The  cumulative  numbers  of  shear  bands  in  Figure  22  indicate 
an  indication  threshold  in  the  vicinity  of  22Z  plastic  strain,  in 
agreement  with  that  predicted  from  the  plastic  instability  theory 
discussed  above  as  well  as  in  Reference  7  and  Appendix  B.  After  the 
initiation  threshold  is  exceeded,  the  sheer  bands  nucleate  gradually  as 
a  function  of  the  increasing  plastic  strain.  The  numbers  appear  to 
fall  into  separate  groups  according  to  tha  explosive  pressure  PCJ. 
indicated  in  the  figure.  The  higher  pressures  also  Induce  a  higher 
strain  rate,  so  the  nudeation  rate  may  be  a  function  of  a  power  of  the 
strain  rate.  The  large  numbers  at  low  strains  for  the  higher  pressure 
tests  may  Indicate  that  most  of  the  Imposed  strain  above  22Z  went  into 
shear  banding  rather  than  into  homogeneous  plastic  flow.  Thus  a 
nucleatlon  process  of  the  following  form  is  indicated 

0.22  (III-25) 

0.22 

to  yield  strength  and  inclusion 


count*.  With  n  •  3  this  expression  for  nuclestlon  rate  agrees  with 

18  31 

that  used  in  our  earlier  work,  ,  ,  for  which  the  corresponding 

expression  for  nuclestlon  in  a  particular  plane,  where  orientation  is 
given  by  the  subscripts  $,6,  is 


Her*  Ecr  is  a  critical  energy,  taken  to  be  the  melt  energy;  dEp/dt  is 
the  rat*  of  increase  of  plastic  strain  energy;  is  the  number  of 
shear  bands  per  unit  volume  with  an  orientation  given  by  the  spherical 
coordinates  $  and  6;  and  ePg  la  the  plastic  shear  strain  in  the  same 
direction.  CQ  is  a  coefficient  with  the  units  of  time  squared,  and 
is  a  factor  representing  the  fraction  of  solid  angle  assigned  to 
the  $6  orientation.  The  factor  dEp/dt  is  proportional  to  tho  product 
of  the  shear  strain  rate  and  the  yield  strength;  hence,  the  nuclaation 
rate  is  proportional  to  the  cub*  of  the  strain  rate  in  SHEAR2,  the 
earlier  version  of  our  shear  banding  model. 

The  shear  band  growth  process  is  studied  by  plotting  R.  as  a 

_p  * 

function  of  e  ,  as  in  Figure  23.  For  the  4340  data,  this  graph  shows 
some  grouping  of  values  for  individual  shots,  but  no  overall  trends. 
The  initiation  slse  may  be  about  *  0.02  cm.  The  total  growth 

produces  about  a  five-fold  increase  in  the  size.  Without  a  clear 
indication  from  the  data,  a  viscous  growth  law  will  be  assumed 

ft  "  CG  H  R 

This  growth  law  preserves  the  exponential  form  of  the  distribution  with 
increasing  strain. 

The  Increase  of  damage  with  plastic  strain  can  also  be  studied 
by  plotting  the  total  shear  band  area  as  a  function  of  the  strain.  The 
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FIGURE  23  SIZE  PARAMETER  R,  FROM  MEASURED  SHEAR  BAND  SIZE  DISTRIBUTIONS 
VERSUS  PLASTIC  STRAIN  BASED  ON  WALL  THICKNESS  MEASUREMENTS  FOR 
SIX  CONTAINED-FRAGMENTING-CYLINDER  EXPERIMENTS  IN  4340  STEEL,  R  40 
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total  area,  a  measure  of  the  damage,  is 


D 


irR2dN 


2trN  R? 
o  1 


(III-27) 


for  tha  exponential  distribution  in  Eq.  (III-18).  The  damage  D,  plotted 
In  Figure  24,  shows  a  clear  trend  toward  Increasing  damage  with  strain, 
and  all  tha  shots  contribute  to  the  same  trend.  The  scatter  at  the 
higher  damage  level  would  probably  be  reduced  by  using  the  imposed 
strain  Instead  of  the  homogeneous  strain. 

The  plots  in  Figures  22,  23,  and  24  represent  the  first  steps 
toward  obtaining  nucleatlon  and  growth  lavs  and  determining  the 
material  constants  for  these  laws.  The  next  steps  are  to  insert  these 
laws  into  SHEAR  3  and  perform  simulations  of  the  experiments  using 
these  constants.  These  simulations  are  under  way. 

2.  Explosive  Pressure  Computation 

In  some  recent  simulations  of  4340  steel  fragmenting  rounds 
for  long  times  (y  100  usee) ,  it  was  found  that  the  test  cylinder 
expanded  beyond  the  dimensions  observed  in  the  experiments.  The 
calculated  pressure  histories  in  these  cases  were  found  to  have  a 
larger  impulse  than  that  expected  from  polytropic  gas  expansion  calcu¬ 
lations.  Therefore,  the  detonation  process  model  was  reexamined.  From 
the  analytical  expression  for  pressure  of  a  running  one-dimensional 
detonation,  we  developed  an  improved  computational  procedure  for  simu¬ 
lating  a  running  detonation.  This  improved  model  is  described  in  the 
following  paragraphs. 

(a)  Analysis  of  One-Dimensional  Detonation.  Behind  the 
detonation  front,  the  explosive  products  (a  polytropic  gas)  may  be 
treated  by  the  method  of  characteristics.  Motion  along  the  C+  charac- 
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FIGURE  24  MEASURED  DAMAGE  VERSUS  PLASTIC  STRAIN  BASED  ON  WALL  THICKNESS 
MEASUREMENTS  FOR  SIX  CONTAINED-FRAGMENTING-CYLINDER  EXPERI¬ 
MENTS  IN  4340  STEEL 


teristlc  is  govsmed  by  Che  following  equation: 


dx 

dt 


u  +  C  - 


(III-28) 


where  x  and  t  are  position  and  time 

u  is  particle  velocity 
C  is  sound  speed 
h  is  Lagrangian  location 

Also  the  particle  velocity  is  defined  as 


u 


h 


and  the  specific  volume  v  is 


v  “  v 

o 


t 


where  v  is  the  initial  specific  volume, 
o 

Replacing  the  definitions  of  u  and  v  in  Eq.  (III-28),  we  obtain 


dh 

dt 


v  C 
o 


(III-29) 


Because  C  and  v  are  constant  along  the  C+  characteristic  in  a  simple 
wave,  Eq.  (HI-29)  can  be  integrated  to  give  the  following  relation 
for  the  rarefaction  fan: 

v  Ct 

h  -  (III-30) 


At  the  front  of  the  detonation,  these  relations  hold: 


yD 


dh 

dt 


Y  +  1 
D 


(III-32) 

(III-33) 


where  Is  the  specific  volume  at  the  C-J  condition 

Y  is  the  polytropic  gas  exponent 

Ccj  is  the  sound  speed  at  the  C-J  condition 

D  is  the  detonation  velocity. 

Now  divide  Eq.  (III-31)  by  Dt  and  replace  v  /D  by  v  ,/C  .  using  Eqs. 

o  cj  Cj 

(III-31)  and  (III-32)  .  Then 


cl  C  h_ 

Ccj  v  Dt 


(III-34) 


Two  additional  facts  from  the  theory  of  polytropic  gases  are  introduced 


Pv' 


P_,v 


Cj  Cj 


C  -  YPv 


(III-35) 

(III-36) 


Combining  these  two  equations,  we  can  obtain 


*-(*) 


Y-l 


(III-37) 


With  this  result  in  Eq.  (III-3I),  we  have 


Vl  _C_ 
V  Ccj 


£2 

v 


2 


h_ 

Dt 


(III-38) 
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Eq.  (III-33)  can  also  bs  ussd  to  relate  the  pressures  to  v/v^y 
Replacing  volumes  with  pressures  in  Eq.  (III-38)  gives 


P 


(III-39) 


Equation  (III-39)  provides  the  pressure  at  any  point  in  the  rarefaction 
fan  behind  the  detonation  front.  For  a  particle  initially  at  a  location 
h,  tbe  detonation  arrives  at  t  -  h/D.  The  pressure  history  for  the 
particle  is  obtained  by  varying  t  from  h/D  on  to  larger  times. 


(b)  Computational  Detonation  Procedure.  For  a  Chapman- Jouguet 
detonation,  the  pressure  should  follow  a  linear  pressure-volume  path 
up  to  the  C-J  point  and  then  expand  down  the  poly tropic  gas  pressure- 
volume  path  given  by  Eq.  (III-35) .  In  addition,  the  C-J  point  should 
be  reached  at  the  time  the  detonation  wave  reaches  the  location.  An 
approximate  procedure  for  meeting  these  requirements  is  outlined  below. 

In  the  finite  difference  calculations,  the  detonation  wave  is 
presumed  to  arrive  when 


t 


(III-40) 


where  Xq  is  the  initial  midpoint  position  of  the  cell 
X^  Is  the  detonation  point. 

Hence  the  detonation  is  forced  along  at  detonation  velocity  through  the 
cells. 

For  each  cell  the  pressure  is  computed  according  to 

P  -  (y  -  l)pE  (III— 4 1 ) 


where  p  is  density  and  E  is  internal  energy. 
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As  chc  density  begins  to  increase  ahead  of  the  detonation  front,  the 
pressure  should  follow  the  linear  path 


P  -  P 


(y  -  V 

cl  vcj  -  vo 


.  pci  (0°  ~  0) 

0  (0ovcj  -  1) 


(III-42) 


To  obtain  this  pressure  in  the  cell,  the  internal  energy  is  prescribed 
by  combining  Eqs.  (111-41)  and  (111-42) 


P 

(Y  -  1)P 


Pc<(0o  -  0) 

Cr  -  i)P!(?0»c:)  -  l) 


When  the  detonation  front  arrives  [Eq.  (111-40)],  the  internal 
energy  is  computed  to  put  the  pressure  on  the  polytropic  expansion 
curve  according  to  Eq.  (111-35) 


P(Y  -  1J 


P(Y  -  1) 


(111-44) 


After  the  polytropic  expansion  curve  is  reached,  the  detonation  process 
is  no  longer  forced.  The  Internal  energy  is  changed  according  to 
the  usual  Pdv  work  relation  and  Eq.  (III-41)  is  used  to  determine  the 
pressure. 


Procedures .  Computations  made  with  each  of  the  foregoing  detonation 
procedures  were  compared.  Pressure-volume  paths  traversed  by  two  cells 
(llth  and  24th)  in  the  finite  difference  computation  are  shown  in 
Figure  25.  The  path  is  initially  along  the  linear  pressure-volume 
Rayleigh  line.  At  the  time  of  detonation,  the  pressure  jumps  up  to  the 
poly tropic  expansion  line  and  then  decays.  Pressure  histories  generated 
at  3.42  cm  (11  cells)  and  7.65  cm  (24  cells)  into  an  explosive  are  shown 
in  Figure  26.  The  analytic  function  shows  an  instantaneous  pressure 
rise  at  detonation  and  a  decay  following  Eq.  (111-39).  By  contrast  the 
finite-difference  computations  show  a  precursor,  a  jump  in  pressure 


FIGURE  26  ANALYTICAL  AND  COMPUTATIONAL  PRESSURE-VOLUME  PATHS  FOLLOWED 
IN  A  DETONATION  PROCESS  IN  LOW  DENSITY  (p  -  1.2  9/cc)  PETN 
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FIGURE  28  COMPARISON  OF  ANALYTICAL  AND  COMPUTATIONAL  PRESSURE  HISTORIES 
IN  LOW  DENSITY  (p  -  1.2  g/cc)  PETN 
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at  tha  time  of  datonatlon  to  a  prassura  balow  P  . ,  and  a  decay  that 

cj 

approaches  tha  analytical  shapa.  Tha  araaa  under  tha  two  curves 
(lapulaa)  are  approximately  equal.  Therefore,  tha  finite-difference 
praasura  history  should  be  sufficiently  accurate  for  determining 
longer  tlaa  phanoaana,  but  cannot  represent  the  wave  front  accurately 
unless  a  such  finer  computational  grid  is  used. 

3.  Analysis  of  RHA  data 

The  shear  band  nuaber  and  size  data  from  the  CFC  experiments 
on  RHA  aaterial  described  in  Section  II  C  have  not  yet  been  reduced, 
and  additional  CFC  experiments  on  this  aaterial  are  under  way.  During 
the  next  year  these  data  will  be  reduced  and  used  to  generate  the  shear 
banding  parameters  in  the  SHEAR3  subroutine. 

C.  Adaptation  of  HEMP  code  for  Use  with  NAG  Material  Models 

The  HEMP  wave  propagation  code  will  be  used  to  simulate  armcr 
penetrations  on  this  project.  In  these  calculations  much  of  the  material 
behavior  will  be  represented  by  the  shear  band  model,  SHEAR  3.  There¬ 
fore,  one  of  our  tasks  is  to  lmert  SHEAR  3  into  HEMP  and  perform  test 
calculations  to  verify  that  the  Insertion  is  correct.  For  the  insertion, 
changes  are  required  in  HEMP,  in  SHZiAR  3,  and  in  the  special  material 
property  data  for  SHEAR  3.  In  cooperation  with  John  French  of  LLL,  we 
have  identified  the  changes  needed. 

The  first  and  most  serious  change  in  HEMP  is  the  Increase  in 
the  number  of  variables  per  cell.  SHEAR  3  requires  6  to  19  variables 

A 

in  addition  to  those  normally  stored  for  e-ich  cell.  Addition  of  more 
variables  requires  Inserting  a  list  of  variables  that  correspond  to  the 
names  used  by  SHEAR  3  and  increasing  the  total  number  of  variables,  MV, 
allowed  per  cell  in  HEMP. 

A 

The  specific  nuaber  depends  on  options  selected  by  the  user. 
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The  CALL  statement  for  SHEAR 3  is  inssrtsd  Into  the  HEMP 
subroutine  celled  VQP  with  ell  the  other  devietor  stress  end  pressure 
celculetions.  VQP  computes  the  strain,  sound  speed,  end  intemel  energy 
in  edditlon  to  the  stresses.  Therefore,  the  second  chenge  required 
for  the  insertion  of  SHEAR  3  is  to  brench  to  the  CALL  ststenent  efter 
the  strein  celculetions.  Stresses  will  be  computed  in  SHEAR  3  end  state¬ 
ments  will  be  edded  in  VQP  for  speciel  celculetions  of  internal  energy 
end  sound  speed  that  ere  consistent  with  the  SHEAR  3  variables.  The 
CALL  statement  to  SHEAR  3  contains  all  the  variables  transferred  between 
SHEAR  3  end  VQP.  Included  ere  the  streeaes,  strein  increments,  the 
variables  describing  sheer  bend  distributions  on  each  plane,  internal 
energy,  yield  strength,  cell  rotation,  end  a  damage  parameter  used  for 
triggering  a  slide  line. 

The  third  change  is  to  provide  for  the  special  input  of 
material  properties  required  by  SHEAR  3.  Currently  in  the  version  of 
HEMP  at  BRL,  the  input  deck  is  praprocesaed  by  the  subroutine  PREHEMP, 
and  a  binary  equivalent  of  this  input  deck  is  written.  Therefore,  to 
provide  for  input  to  SHEAR  3,  we  must  modify  the  read  statements  in 
SHEAR  3  to  accept  this  binary  input.  The  data  for  SHEAR  3  will  follow 
the  normal  HEMP  input.  The  density,  the  shear  modulus,  and  the  bulk, 
modulus  will  be  Included  in  the  normal  HEMP  material  input,  and  only 
special  shear  banding  variables  will  be  read  by  SHEAR  3. 

The  changes  in  SHEAR  3  required  by  the  Insertion  also  include 
rewriting  the  read  statements  for  the  material  properties.  To  change 
SHEAR  3  from  the  dyn-sec-cm-g  system  of  units  to  the  Mbar-usec-css- g 
system  used  in  HEMP  should  require  little  recoding. 

Finally,  the  special  material  property  data  for  SHEAR  3  are 
being  transformed  to  the  HQCP  system  of  units. 

The  above  changes  are  currertly  under  way.  Furthermore,  the 
changes  required  to  insert  other  NAG  material  models  such  as  BFRACT 
(dynamic  brittle  feature),  DFRACT  (dynamic  ductile  fracture),  and 
DFRACTS  (quasi-static  ductile  fracture)  will  be  very  similar  to  those 
required  to  Insert  SHEAR3.  The  SHEAR3  insertion  is  therefore  being 
accomplished  In  a  manner  designed  to  make  later  insertion  cf  the  other  sub¬ 
routines  as  efficient  as  possible. 
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0.  Compute  tlonal  Simulations 


1.  Background  and  Scop* 

Tha  ultimata  goal  for  this  program  la  to  ba  abla  to  computa¬ 
tionally  slmulata  apaced  armor  panatration  with  long  rods,  and  to  ba 
abla  to  pradict  the  fragmant  environment  behind  tha  last  plate.  To 
attain  this  goal  will  require  combining  adequate  micros true tural  material 
failure  models  with  state-of-the-art  computer  codes.  It  is  therefore 
important  aarly  in  the  program  to  assess  the  feasibility  of  our  approach. 

The  simplest  possible  case  (but  one  rarely  seen  in  practice 

with  RHA)  is  normal  impact  with  plugging  of  the  target,  such  that  the 

fragment  environment  behind  the  plate  consists  of  a  single  or  fragmented 

plug  plus  penetrator  fragments.  During  this  first  year's  effort  we  have 

9  10 

performed  computational  simulations  of  this  case  ,  and  have  confirmed 
that  it,  at  least,  is  amenable  to  the  combined  microstruc tural  failure 
model/computational  approach.  In  the  following  paragraphs  we  describe 
some  of  these  calculations. 

To  calculate  the  plugging  process,  several  factors  must  be 
modeled  correctly.  The  first  phase  in  the  formation  of  the  plug  is  the 
shear-band  nudeatlon  and  growth.  The  second  phase  is  the  coalescence 
of  the  shear  bands  into  a  macroslip  surface.  Here  it  is  necessary  to 
model  both  the  threshold  conditions  at  which  this  coalescence  takes 
place  and  the  location  of  the  coalesced  surface.  The  third  phase  is 
the  motion  of  the  plug  as  it  is  pushed  out  of  the  target. 

In  these  calculations  the  older  SHEAR2  model  was  used  to 
calculate  the  microdamage,  since  SHEAR3  was  still  under  development. 

2.  Unsipping-Slide-Line  Logic 

The  purpose  of  the  unzipping  slide  line  is  to  model  the 

coalescence  of  microshear  bands  into  a  macroshear  band  or  slide  line. 

In  the  current  calculations  the  slide  line  lies  along  a  predetermined 

position  but  eventually  the  position  of  the  slide  line  must  also  be  de- 

q 

termined  by  the  mlcroshear-band  model.  The  strategy  of  the  unzipping 
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process  is  co  allow  chs  microdamage  to  inersaas  until  a  threshold  level 
of  damage  is  reached.  When  any  of  the  computational  cells  around  a 
nodal  point  lying  on  the  slide  line  meet  this  criterion,  the  nodal 
point  changes  from  a  standard  nodal  point  to  a  slide-line  nodal  point. 

In  a  standard  calculation,  a  nodal  point  that  is  not  a 
.  boundary  point  is  surrounded  by  four  cells;  A  in  Figure  27  is  such  a 

point.  The  stress  in  these  four  cells  will  contribute  to  the  force  acting 
on  the  nodal  point.  This  force  is  then  used  to  calculate  the  new 
accelerations,  the  new  velocities,  and  the  new  coordinate  positions. 

A  slide-line  calculation  treats  nodal  points  differently 
frost  the  standard  calculation.  Here  we  have  two  distinct  ways  of  cal¬ 
culating  the  forces  acting  on  a  nodal  point,  depending  on  whether 
the  nodal  point  is  a  master  or  slave. 

A  series  of  master  nodal  points  will  determine  the  position  of 
the  boundary,  and  the  slave  nodal  points  will  slide  along  that  boundary. 

For  a  master  nodal  point  that  is  not  a  boundary  point,  the  stresses  of 
three  cells  will  contribute  to  the  force  acting  on  the  nodal  point  and 
hence  determine  its  motion.  On  the  other  hand,  the  slave  nodal  point 
derives  its  motion  from  stresses  in  only  the  neighboring  two  cells; 

ithe  position  of  the  nodal  point  is  then  adjusted  to  slide  along  the 

master-defined  slide  line.  In  Figure  27,  B  represents  a  master  r.odal 
point  and  C  represents  a  slave  nodal  point.  Special  consideration  must 
be  given  if  any  of  these  points  are  either  the  first  or  last  points  of 
the  slide  line. 

In  a  calculation  in  which  the  slide  line  is  unzipped,  a  nodal 
point  first  must  be  treated  as  a  regular  nodal  point — that  is,  as  A 
in  Figure  27 .  If  any  computed  value  in  the  neighboring  cells ,  such  as 
plastic  strain  or  present  shear-band  damage,  meets  a  threshold  condition,  j 

the  single  nodal  point  then  separates  into  two  separate  points  and  the 
forces  are  calculated  as  in  B  and  C  of  Figure  27 . 

•  At  any  time  in  the  calculation,  the  slide  line  will  consist  1 

a 

of  a  combination  of  zipped  and  unzipped  nodal  points.  As  the  calculation 
progresses,  more  nodal  points  are  unzipped  until  a  complete  slide  line 
is  formed. 
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3.  Computational  Setup 


A  series  of  computations  were  performed  for  a  normal  rod  impact, 
with  various  treatments  of  the  slide  line.  Figure  28  shows  the  computa¬ 
tional  layout  for  a  1.27-cm-thick  4340  steel  target  impacted  by  a  4340 
steel  rod  1.03  cm  in  diameter  and  2.06  cm  long.  In  all  the  calculations, 
the  impact  velocity  was  750  m/sec,  which  experiments  indicated  was  near 
the  ballistic  limit. 

Figure  29  shows  a  cross  section  of  the  target  from  such  an 
experiment  f impact  velocity  *  790  m/sec).  It  is  clear  that  the  shear  bands 
forming  the  plug  occur  near  the  rod  edges  at  impact,  thus  making  the 
assumed  position  of  the  incipient  slide  line  shown  in  Figure  28  appropriate. 

In  these  calculations,  we  computed  shear  strains  in  six  orientations. 
Three  of  the  orientations  are  for  band  planes  normal  to  the  following  directions : 
(1)  axial,  (2)  radial,  and  (3)  circumferential.  The  other  three  orientations 
are  for  band  planes  that  are  at  45°  between  the  following  directions:  (4)  axial 
and  radial,  (5)  axial  and  circumferential,  and  (6)  circumferential  and 
radial.  Only  the  first,  second,  fourth,  and  fifth  orientations  were  ex¬ 
pected  to  be  important  in  this  case;  therefore,  we  calculated  the  shear- 
band  distributions  only  in  these  directions,  but  continued  to  calculate  the 
directional  plastic  strain  in  all  six  orientations.  In  this  calculation, 
band-size  distributions  are' represented  by  a  series  of  discrete  size  groups. 
(The  model  option  with  an  analytical  form  for  the  size  distribution  was  not 
used.)  For  each  of  the  orientations,  there  are  six  size  groups. 

The  purpose  of  these  calculations  was  to  examine  the  unzipping 

sliJe  line  logic.  Since  the  SNAG  parameters  for  4340  steel  are  still  being 

determined,  we  used  instead  the  parameters  previously  determined  for  a 

33 

high  fragmentation  steel.  For  this  material,  the  maximum  bind  size  at 
nucleafion  is  0.07  cm.  The  radii  of  the  bands  are  initialized  so  that  the 
ratio  of  successive  intervals  between  radii  is  equal  to  1.4.  Hence,  if 
the  maximum  radius  is  0.07  cm,  the  next  size  radius  group  is  0.05  cm,  and 
so  on.  The  threshold  equivalent  plastic  shear  strain  for  shear  band 
nucleation  is  C..7. 


The  values  of  all  the  parameters  used  are  summarized  in  Table 
4.  In  addition  to  the  usual  variables  carried  per  cell,  63  additional 
variables  were  used.  ^ 
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FIGURE  28  INITIAL  CONDITIONS  FOR  CALCULATION  OF  A  STEEL  PROJECTILE 
OF  VELOCITY  0.75  km/sec  IMPACTING  A  STEEL  TARGET 

The  slide  line  is  defined  by  10  nodal  points  and  Is  unzipped  during  the  * 

calculation,  r 
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Table  4 


SHEAR2  MATERIAL  PROPERTIES  USED  IN  COMPUTATIONAL 
SIMULATIONS  OF  4340  STEEL  ROD  IMPACTING  NORMALLY  A  4340  STEEL  TARGET 


Ruelaaclon  equivalent  plastic,  threshold  strain, 


Maxima  band  size  at  nudeatlon,  R 


Ratio  of  successive  Intervals  between  band  sizes 


Nucleatlon  size  paramter,  Rj,  froa  Eq.  (III-18) 


Nudeatlon  rate  paramter,  A^,  froa  Eq.  (II I- 25) 


Growth  rate  paramter  froa  Eq.  (III-26) ,  C 


Ratio  of  band  jog  to  band  radius,  b 


Number  of  orientations  used,  NANG 


Number  of  size  groups  assigned  to  each  orientation. 


Orientation 


Number 


0.07  ca 


0.011  cm 


2  -3 

0.001  sec  cm 


Critical  internal  energy  from  Eq.  (III-251),  Ecr  7  x  10  erg/g 


N 
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As  shown  in  Figure  28,  the  slide  line  was  placed  in  the  target 
directly  under  the  edge  of  the  penetrator  to  form  the  boundary  of  the 
plug.  Whether  this  slide  line  is  unzipped  as  the  calculation  progresses 
depends  on  whether  the  specified  threshold  condition  has  been  reached. 

In  this  calculation  the  slide  line  is  defined  by  10  nodal  points. 

To  examine  the  unzipping  slide  line  behavior,  we  used  two 
threshold  criteria  for  unzipping.  The  first  criterion  was  a  critical 
value  of  the  equivalent  plastic  strain.  By  systematically  varying  this 
value,  we  produced  various  amounts  of  unzipping  and  tested  the  stability 
of  the  unzipping  logic.  These  calculations  are  described  in  Reference  10. 

However,  this  unzipping  criterion  is  unrealistic  because  it  is  not 
directly  connected  to  the  amount  of  microscopic  shear  band  damage. 

Therefore,  after  the  unzipping  logic  had  been  tested,  a  second  criterion 
was  used,  namely  a  critical  value  of  the  damage  parameter  x,  where  x  is 
the  fraction  of  a  material  cell’s  strength  that  has  been  removed  by  shear 
banding.  The  results  of  some  of  these  calculations  are  discussed  in 
the  following  paragraphs. 

4.  Selected  Computational  Results 

Two  computations  are  described  here.  The  first  deals  with  an 
unzipping  slide  line  for  which  the  unzipping  is  governed  by  a  critical 
plastic  strain.  The  second  calculation  deals  with  a  slide  line  that 
uses  a  damage  parameter  calculated  from  the  shear-banding  subroutine  to 
trigger  the  unzipping  process.  The  results  of  this  calculation  are 
compared  with  the  results  of  the  impact  shown  in  Figure  29. 

(a)  Plastic  Strain  Criterion  for  Slide  Line  Activation.  The 
results  of  a  calculation  in  which  a  threshold  plastic  strain  of  10%  was 
used  to  trigger  the  unzipping  of  a  slide  line  are  shown  in  Figure  30. 

The  process  is  shown  in  three  successive  stages  of  development.  Here 
the  slide  line  was  defined  by  10  nodal  points;  no  nodal  point  slid  until 
one  of  its  neighboring  cells  reached  the  10%  critical  plastic  strain 
threshold.  The  unzipping  of  the  slide  line  (complete  formation  of  the 
plug)  took  about  10  ysec  or  120  computational  cycles. 
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This  calculation  was  one  of  the  series  reoorted  in  Reference  10 
to  test  the  unzipping  logic,  and  is  reproduced  here  to  show  several 
features.  First,  when  the  penetrator  overmatches  the  target,  the  pro¬ 
jectile  nose  does  not  mushroom  enough  to  impact  the  target  face,  and 
only  one  active  slide  line  is  required.  Second,  the  macroshear  band 
advances  in  a  reasonable  way  through  the  target  and  releases  the  plug 
without  computational  difficulty.  Third,  the.  behind-the-armor  fragment 
environment  would  be  described  completely  by  the  unfragmented  portion  of 
the  plug  plus  the  shear  band-produced  fragmentation  of  the  plug  and 
penetrator  nose. 

(b)  Shear  Band  Damage  Criterion  for  Slide  Line  Activation. 
Following  the  above  debugging  calculations  to  examine  the  unzipping 
slide  line  behavior,  we  performed  a  simulation  in  which  the  unzipping 
criterion  was  a  critical  value  of  the  parameter  t 

t  -  f  it  2  N^3, 

which  is  related  to  the  percentage  of  microscopic  shear  band  damage  in 

the  cell.  N  is  the  number  of  bands  per  unit  volume,  R  is  the  band 

radius,  and  the  subscript  i  refers  to  the  band  size  group.  The  critical 

* 

value  of  t  chosen  for  the  calculation  was  2.0.  The  resulting  deformed 
grid  at  7.98  usee  after  impact  is  shown  in  Figure  31(a).  The  calcula¬ 
tion  was  stopped  at  7.98  ysec  when  the  projectile  bulge  began  to  impact 
the  target  shoulder,  thus  necessitating  a  second  slide  line  between 
the  front  of  the  target  and  the  side  of.  the  projectile,  as  discussed 
above.  At  this  time,  approximately  the  front  third  of  the  rod  was 
calculated  to  have  been  fragmented  by  shear  banding,  the  rod  still  had 
44Z  of  its  original  momentum,  and  the  velocities  of  the  front  and  rear 
of  the  rod  had  dropped  to  88  ra/sec  and  607  m/sec,  respectively.  The 


The  selection  of  this  value  was  an  error  that  was  undetected  until 
after  the  calculation  was  performed.  The  desired  value  was  1, 
corresponding  to  complete  local  fragmentation  due  to  shear  banding.  The 
effact  of  this  larger  value  of  t  will  be  to  delay  the  unzipping  somewhat. 
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(•)  CALCULATION  7.98  psec  AFTER  IMPACT 


j  PROJECTILE 


mm  Fragmented 
“““  Region 

(b)  COMPARISON  OF  COMPUTED  WITH  OBSERVED  PENETRATION 
BEHAVIOR 

In  the  calculation  the  projectile  still  hat  44%  of  its  original  momentum, 
but  the  front  portion  it  fragmented  by  shear  banding.  The  calculation 
was  stopped  because  the  shoulder  of  the  projectile  has  struck  the  target 
face,  thus  requiring  another  slide  line. 

MP-7893-72 

FIGURE  31  UNZIPPING  SLIDE-LINE  CALCULATION  WITH  THRESHOLD  SHEAR  BAND 
DAMAGE  PARAMETER,  r  -  2.  AND  COMPARISON  WITH  EXPERIMENT 
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central  deflection  of  the  target  in  0.22  cm  and  the  protrusion  at  the 
back  Is  0.03  cm.  The  slide  line  has  only  been  activated  for  a  short 
distance  into  the  target  as  shown.  The  diameter  of  the  slide  line 
surface  has  expanded  from  the  original  1.03  cm  to  1.28  cm. 

These  computational  results  were  compared  with  the  experimental 
results  of  Shot  6802-4  shown  in  Figure  29  to  gain  insight  into  the  pene¬ 
tration  process;  the  comparison  is  shown  in  Figure  31(b).  The  comparison 
can  only  be  quantitative  because  the  calculation  was  stopped  before  the 
projectile  was  completely  fragmented  or  had  come  to  rest,  the  shear  band 
parameters  were  for  a  different  steel  than  4340,  and  an  unrealistically 
large  value  of  T  was  inadvertently  used.  Nevertheless,  we  consider  the 
comparison  very  promising.  The  shear  banding  and  formation  of  a  macro- 
band  or  plug  appears  to  be  occurring  correctly  in  the  calculation,  and 
the  front  surface  crater  and  rear  surface  bulge  are  developing  as 
observed  experimentally.  The  mushrooming  of  the  rod  nose  and  impact 
with  the  target  face  also  seem  to  occur  both  computationally  and  exper¬ 
imentally  at  this  impact  velocity.  The  comparison  shows  that  a  double 
slide  line  is  required  to  correctly  represent  the  penetration  in  some 
cases,  and  will  certainly  be  required  in  most  oblique  impact  cases. 

During  the  next  year  we  plan  to  perform  more  such  calculations 
using  the  SHEAR3  subroutine  in  place  of  SHEAR2;  both  4340/4340  and 
DU/RHA  normal  impact  experiments  over  a  range  of  impact  velocities  will 
be  simulated  to  test  the  model's  ability  to  reproduce  the  main  features 
of  normal  penetration  and  fragmentation.  In  addition,  work  will  begin 
on  simulation  of  oblique  impacts. 


IV  CONCLUSIONS  AND  FUTURE,  PLANS 


The  key  conclusion  Co  dace  Is  chat  adiabatic  shear  b,  nding  is  very 
probably  Che  dominant  mechanism  for  plugging  and  fragmentation  of  RHA 
and  for  erosion  of  the  DU  penetrator  nose.  Furthermore,  the  computational 
shear  band  model  SHEAR3  should  provide  adequate  modelling  of  this 
behavior.  We  will  test  these  conclusions  further  during  the  next  year 
by  performing  higher  impact  velocity  phenomenology  tests  and  additional 
computational  simulations. 


In  the  upcoming  phenomenology  experiments  and  computations,  we  will 
concentrate  on  the  effects  of  impact  obliquity  and  yaw,  since  most  of 
the  quantltatlveresults  to  date  concern  only  normal  impacts. 

Tensile  failure  by  ductile  void  nucleation,  growth,  and  coalescence 
is  probably  the  second  key  failure  mode,  governing  spall  and  petalling 
of  the  target  and  fracture  of  the  rear  portion  of  the  DU  rod.  However, 
further  phenomenology  experiments  are  required  to  confirm  this  conclusion. 
NAG  models  have  been  developed  and  are  ready  for  application  to  the 
phenomenology  experiment  results. 


We  conclude  further  that  the  program  organization  and  milestones 
diagrammed  in  Figure  1  are  still  valid  and  feasible,  and  we  shall 
continue  to  follow  that  approach  during  the  next  year. 


Specific  tasks  for  th  e  second  year  of  the  contract  were  determined 
during  a  meeting  at  AMMRC  in  August  1979.  Depleted  uranium  (DU), 
tungsten,  and,  possibly,  electro-slag-remelted  (ESR)  steel  are  to  be 
procured.  Targets  and  penetrators  from  experiments  performed  at  BRL 
and  AMMRC  are  to  be  obtained  and  examined  to  determine  the  dominant 
fracture  modes.  Additional  phenomenology  experiments  with  DU  rods  and 
RHA  targets  will  be  conducted  at  SRI  and  at  BRL  with  appropriate  obliquity, 
velocity,  and  yaw  to  further  define  the  fracture  modes  of  both  elements. 
Contained-fragmenting-cylindei  tests  and  other  fracture-characterizing 
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experiments  will  be  conducted  by  SRI  on  RHA  end  DU.  (Desciaions 
concerning  specific  experiments  will  be  based  on  the  results  of  the 
phenomenology  tests.)  The  improved  shear-band  model,  SHEAR2,  will  be 
applied  to  the  4340  steel  data,  tested  by  simulation  of  a  penetration 
experiment  into  a  4340  steel  target,  and  then  applied  to  the  RHA  and 
DU  data  to  represent  plugging  and  penetrator  nose  erosion.  SHEAR3  will 
be  implemented  in  the  HEMP  code  arid  will  be  exercised  with  double  slide 
lines.  These  planned  tasks  are  discussed  in  more  detail  below. 

A.  Materials  Selection  and  Phenomenology  Experiments 

BRL,  AMMRC,  and  SRI  have  jointly  decided  to  study  a  uranium-3/ 4 

ee  A 

vtZ  titanium  alloy  first.  The  uranium  used  is  depleted  uranium,  * 

the  fissionable  isotope  having  been  removed. 

During  the  past  few  months,  a  supply  of  the  alloy  was  selected, 
the  material  was  fabricated,  and  the  necessary  specimens  were  machined 
and  shipped  to  SRI.  The  specimens  were  machined  only  from  rods 
containing  centerline  porosity  of  less  than  40  jim  in  diameter. 

1,  DU  Rods  Penetration  Phenomenology  and  Characterization 

Experiments 

Ten  4.7-g  hemispherical-nosed  rods,  0.125  inches  in  diameter 
and  1.250  inches  long  (L/D  •  10),  were  machined  from  the  bar  stock. 

The  SRI  ballistic  test  facility  wij 1  be  used  to  fire  these  rods  at 
previously  characterized  Rc22  RHA  targets  at  oblique  angles  to  determine 
the  failure  phenomenology  of  U-3/4  wt?  Ti  alloy.  In  addition,  15  65-g 
hemispherical-nosed  rods,  0.3023  Inches  in  diameter  and  3.023  inches 
long  (L/D  »  10),  were  machined  from  the  same  1.4-inch-diameter  bar 
stock.  These  rods  will  also  be  fired  at  oblique  R222  RHA  targets, 
using  our  light  gas  gun,  to  study  phenomena  scability.  After  the 
phenomenology  experiments  have  determined  what  failure  or  deformation 
modes  require  quantitative  characterization,  thesa  rods  may  also  be 
used  for  DU  characterization  experiments. 

Both  the  4.7-g  and  the  65-g  rods  were  machined  from  the  bar 
stock  in  such  a  manner  as  to  avoid  measurable  centerline  porosity. 
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In  addition  to  thaaa  phenomenology  experiments  to  be  performed 
at  SRI,  additional  DU/RHA  experiments  will  be  performed  et  higher  Impact 
velocities  at  BRL  under  the  direction  of  W.  Gillich.  The  purpose  of 
these  tests  will  be  to  obtain  behind-the-armor  fragmentation  data  as 
well  as  phenomenology  information. 

2.  DU  Tubes  for  Confined-Fragmenting-Cylinder  Experiments 

If  adiabatic  shear  banding  is  confirmed  to  be  the  dominant 
failure  mechanism  in  U-3/4  wtX  T1  during  ballistic  impact,  a  series  of 
confined-fragmenting-cylinder  experiments  will  be  conducted  to  determine 
the  necessary  parameters  for  the  SKEAR3  computational  model.  A  total 
of  12  tubes,  1.280  inches  outside  diameter,  0.880  +  0.0005  inch  inside 
diameter,  and  4.0  inches  long,  have  been  nnchined  from  the  same 
1.4- inch-diameter  bar  stock.  These  tubes  were  difficult  to  machine 
because  of  the  nature  of  the  alloy  and  the  dimensional  tolerances 
required . 

All  12  tubes,  all  15  65-g  projectiles,  and  all  10  4.7-g 
projectiles  have  been  received  by  SRI. 

3.  Tungsten  Alloy  Penetrators 

Discussions  are  now  under  way  with  BRL  and  AMMRC  personnel  to 
arrange  for  SRI  to  receive,  for  metallographic  examination,  penetrator 
and  rod  fragments  from  a  series  of  penetration  experiments  planned  at 
BRL  and  AMMRC.  These  experiments  will  be  discussed  more  thoroughly 
in  future  reports. 

B.  Material  Characterization  Experiments 

1 .  Contained-Fragmenting-Cylinder  (CFC)  Experiments 

During  the  next  year  several  additional  CFC  experiments  will 
be  performed  on  RHA  material,  and  experiments  may  be  initiated  with 
DU  material.  The  SHEAR3  model,  as  modified  by  correlation  with  the 
4340  steel  data,  will  then  be  applied  to  both  the  RHA  and  DU  data. 
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2.  Symcetric  Taylor  Tests 


Several  symmetric  Taylor  teats  on  RHA  and  DU  rods  will  be 
performed,  and  the  data  will  be  analyzed  with  the  help  of  2-D  code 
calculations  to  determine  the  dynamic  yield  functions. 

3.  Smooth  Bar  Tensile  Tests 

We  will  initially  attempt  to  obtain  smooth  bar  tensile  test 
data  for  DU  from  the  literature.  If  such  data  are  unavailable  or 
ambiguous,  we  may  perform  several  such  tests  at  SRI. 

4.  Plate  Impact  Experiments 

Plate  impact  experiments  to  obtain  the  NAG  parameters  for 
dynamic  ductile  void  nucleation,  growth,  and  coalescence  in  (composition) 
tungsten  alloy  have  been  previously  performed  by  Gerald  Moss  at  BRL, 
are  currently  under  way  at  BRL  for  our  baseline  RHA  material,  and  are 
planned  during  the  upcoming  year  for  DU  material.  The  DU  material  is 
from  the  same  batch  described  above  for  the  rod  and  cylinder  tests. 

C.  Computational  Model  Adaptation 

1.  SHEAR3 

As  discussed  above,  the  SHEAR3  subroutine  for  adiabatic  shear 
banding  is  currently  being  debugged  and  used  to  correlate  the  results 
of  computational  simulations  with  4340  steel  CFC  data.  No  further 
development  of  SHEAR3  is  foreseen  during  the  next  year;  instead,  we 
will  concentrate  on  documentation  and  application  of  SHEAR3  to  DU/RHA 
penetration  problems. 

2.  DFRACT  and  DFRACTS 

DFRACT,  the  NAG  nubroutine  for  dynamic  ductile  void  growth, 
will  be  used  in  its  current  form  to  correlate  BRL's  plate  impact 
experimental  data  for  Tungsten,  RHA,  and  DU,  so  as  to  obtain  the  dynamic 
NAG  material  parameters  for  these  materials.  No  additional  code 
development  of  DFRACT  is  foreseen. 


DFRACTS,  the  NAG  subroutine  for  quasi-static  ductile  void 
growth,  has  yet  to  be  Inserted  in  wave  propagation  codes  and  debugged. 
These  steps  will  be  delayed  until  the  phenomenology  experiments  have 
established  the  importance  of  this  damage  mode  in  armor  petalling  or 
penetrator  fracture. 

D.  Generation  and  Testing  of  Model  Parameters 

As  discussed  above,  micro structural  failure  data  will  be  available 
during  the  next  year  for  4340  steel,  RHA  material,  DU,  and  tungsten. 
These  data  will  include  results  from  quasi-static  tensile  tests,  plate 
impact  tests,  CFC  tests,  and  symmetric  rod  impact  tests.  These  data 
will  be  correlated  with  computational  simulations  to  generate  and  test 
model  parameters  for  DFRACT,  SHEAR3,  and,  possibly  DFRACTS.  During  the 
next  year  we  plan  to  complete  this  program  for  SHEAR3  using  4340  steel 
and  RHA  material,  whereas  the  remaining  work  would  be  completed  during 
the  third  year. 

E.  Adaption  of  HEMP  Code  for  Use  with  NAG  Material  Models 

During  the  next  year  the  subroutine  SHEAR3  will  be  inserted  in 
HEMP,  debugged,  and  run  on  penetration  problems  with  double  slide  lines 
and,  if  necessary,  rezoning.  We  will  begin  insertion  of  the  subroutines 
DFRACT  and,  if  necessary,  DFRACTS.  These  insertions  will  be  documented 
in  cooperation  with  BRL  and  AMMRC  personnel. 

F.  Computational  Simulations 

Successfully  implementing  SHEAR3  in  HEMP  and  testing  with  double 
slide  lines  and  rezoning  will  occupy  a  major  portion  of  the  second 
year's  effort.  The  computational  simulations  involved  in  this  effort 
will  be  of  two  types.  In  the  first,  computational  simulations  will 
be  attempted  of  normal  impacts  against  single  RHA  plates  ranging  from 
near  the  ballistic  limit  to  cases  where  the  target  plate  is  significantly 
overmatched.  We  will  compare  the  results  to  phenomenology  experimental 
data,  and  study  the  ability  of  SHEAR3  to  successfully  predict  penetrator 
and  target  fragmentation .and  the  behind-the-armor  fragment  environment. 


85 


la  the  second  computational  simulation  effort,  we  will  begin 
simulating  cases  of  varying  impact  obliquity  and  yaw.  In  these  cases, 
the  two-dimensional  nature  of  the  codes  will  restrict  us  to  plate-on- 
plate  simulation.  Previous  work  at  LLLJ  and  BRL  has  shown  that  the 
mein  features  of  oblique  rod  impact  can  be  simulated  with  the  plate- 
ott-plate  approximation. 
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APPENDIX  A 

THE  SHEAR 3  SUBROUTINE  FOR  ADIABATIC  SHEAR  BANDING 
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1.  Modification!  to  SHEAR2 


Th«  shear  band  nodal  haa  baan  axtansivaly  nodlflad  slnca  tha  docu- 
18  31 

mentation  of  SHEAR2.  *  Tha  ahaar  band  alaa  distributions  wars  changed 
fro*  a  sarias  of  dlscrata  points  to  an  analytical  fora  (Rtf.  2),  and  tha 
stress  relaxation  associated  with  dasaga  was  coapletaly  raderived  (Raf.  4). 
Additional  Modifications  ware  made  in  recant  aonths  to  incorporate  tha  new 
mideation  and  growth  relations  and  to  reduce  tha  nunber  of  special 
variables  par  call  required  by  tha  subroutine.  This  reduction  was 
accomplished  by: 

a  Eliminating  XH3,  BN,  and  FLO)  (which  represent  an  indicator, 
tha  total  number  of  bands  in  all  orientations  in  each  call, 
and  Che  plastic  work). 

•  Providing  an  option  for  storage  of  plastic  strains  in  each 
orientation.  The  strains  are  not  required  for  the  calculation 
but  .ire  of  great  interest  in  assessing  the  relative  potential 
for  damage  of  each  orientation.  This  option  permits  the  uier 
to  avoid  storing  these  strains. 

»  Reducing  the  shear  band  directions  from  6  to  4  on  planes: 

(1)  normal  to  x  or  y,  (2)  at  45*  between  x  and  y,  (3)  at  45* 
between  y  and  a,  and  (4)  at  45*  between  3  and  i.  This  change 
actually  does  not  reduce  the  information  content  from  that  in 
the  six  orientations  previously  used.  Vith  the  current  nudes- 
tion  and  growth  process,  the  bands  normal  to  x  and  to  y  are 
always  identical,  so  one  distribution  can  represent  both 
orlentutlons .  There  are  no  bands  normal  to  the  s  direction 
(circumferential  in  axlsyametric  problems)  so  that  orientation 
is  unneeded. 

•  Providing  an  option  for  user  selection  of  shear  band  orientations. 
The  user  can  choose  to  not  activate  any  of  the  orientations. 

With  the  present  configuration  of  SHEAR3,  the  additional  variables 
required  beyond  those  for  a  standard  elastic-plastic  model  are  as  follows: 

TAU,  the  damage  parameter,  to  be  used  for  triggering  slide  lines 

EP,  EP 
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ROT,  the  call  rotation 

tha  plaatle  strain  In  aach  orlantatlon  (0  to  4  variables) 

threa  variables  (N  ,  R. ,  and  R _ )  for  aach  orientation  activated. 

O  1  mmX 

Tha  call  statements  for  SHEAR3  ara  similar  to  thosa  for  SHEAR2.  Tha 
first  call  Is  used  to  raad  In  data  and  initialise  tha  aubroutlna:  this 
call  la  from  LATOUTT  in  TROTT.  Only  tha  parameters  NCALL  and  IN  are 
required  for  this  Initialising  call.  Tha  second  call  stateamnt  is  in 
SWEEPT  of  TROTT.  This  call  is  for  tha  stress  and  damage  calculations 
and  requires  all  tha  call  variables.  Note  that  all  data  passes  through 
the  call  statement:  there  are  no  COMMON  variables.  The  call  and  associated 
statements  in  SWEEPT  are: 

492  SXXW-SXX(LM) 

STTW-STY(LM) 

TZTV-TZT(LM) 

Ftf-P(LM) 

LS  -  2 

IP  (N0D( IN, INPRINT)  .  EQ.  0)  LS  -  3 

CALL  SHEARS  (LS,5»MAT,E,J,8XXW,SYTV,TXYH,FW,C0M(LM+20) , 

1  0W,0(LM)  ,DT,EU,E(LM)  ,C0M(LM+21)  ,EXXH,ETTH, 

2  EXTH,F,TT(LM)  ,TH(U() , -ALFA,  ISC, COM(LM+22)) 

IF  (C0MOM4-20)  .CT.  0.)  IB(LM)-3 

szzh— sexh-sttv 

TH(LM)-TH(LM)-ALFA 
00  TO  600 

The  first  four  statements  bring  the  deviator  stresses  and  pressure  front 
the  akaln  arrays  in  preparation  for  the  call.  LS  la  set  to  2  for  normal 
stress  calculations  and  to  3  for  calculations  and  printing..  Following 
the  call,  the  Indicator  IH  Is  set  to  3  if  shear  banding  has  occurred: 

GQM(LK  <6  20)  la  the  main  storage  for  T  ■  TAU,  a  damage  parameter.  The 
rotation  factor  TH  la  Incremented  by  ALFA,  the  differential  rotation, 
which  Is  positive  clockwise  in  TROTT.  TH  is  positive  counter  clockwise. 

The  stresses  and  strains  are  positive  in  tension;  the  pressure  is  positive 
in  compression.  The  variables  in  the  call  statement  are  listed  in  Table 
A-l. 

5H1AR3  has  been  used  with  the  dyn-cm-g-sec  system  of  units  exclusively. 
He  have  searched  through  the  subroutine,  hunting  for  constants  that  would 
require  changes  if  different  systems  are  used.  For  the  Mbar-cm-g-psec 
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TABLE  A-l 


VARIABLES  IN  THE  CALL  TO  SHEAR  3 


Variably 

NCALL 


IH 

M 

K 

J 

SI 

ST 

SZ 

SXT 

P 

TAU 

DH 

BOLD 

DTO 

EH 

SOLD 


_ Meaning _ 

Indicator  for  tho  type  of  calculations  raquirad 
in  SHEAR  3 

0,1  Raad  data  for  ona  natarlal  and  initialise 

2  Calculate  only 

3  Calculate  and  print 

4  Print  only 

Data  input  file  number  (usually  5) 

Material  number 

Call  number  in  the  K  direction 
Call  number  in  the  J  direction 

2 

Deviator  stress  in  the  X  direction  (dyn/cn  ) 

2 

Deviator  stress  in  the  T  direction  (dyn/cm  ) 

2 

Deviator  stress  in  the  Z  direction  (dyn/cm  ) 

2 

Shear  stress  (dyn/cm  ) 

2 

Pressure  calculated  for  the  cell  (dyn/cm  ) 

Damage  (E  NR^)  in  the  cell  calculated  on  the  last 
cycle;  ranges  from  0  to  1  (See  BFft(M,32) 

New  incoming  density  (g/cm  } 

Cell  density  determined  on  the  last  cycle  (g/cm^) 
Tima  step  (sec) 

Internal  energy  of  the  cell  in  the  present  cycle 

(•rg/g) 

Energy  of  cell  in  previous  cycle  (erg/g) 
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TABU  A-l 


Variable 

KP 

n 

IT 

KZ 

m 

p 

THL 

IDT 

DBOT 

ESC 

ESC(M, 1) 

ESC(M,2) 

ESC(M,3), 
ESC(M*4  i 

ESC(M,5) 

ESC(M,9) 

CM 


(Continued) 


M»inln| 

Total  offoctlvo  plastic  strain 

Strain  increment  In  tha  Z  direction 

Strain  increment  In  tha  T  direction 

Strain  Increment  In  the  Z  direction 

Tensor  shear  strain  lncrsasnt 

Thermal  softening  factor;  varies  frosi  0  to  1 

Yield  value  calculated  on  the  last  cycle  [remains 
unchanged  (1)  If  the  yield  has  not  been  exceeded* 
or  (2)  If  it  has  been  exceeded  and  work  hardening 
Is  saro]  (dyn/cm2) 

Total  angle  of  cell  rotation  as  of  the  previous 
cycle 

Additional  angle  of  cell  rotation  since  the  last 
cycle 

Equation  of  state  variables  [array  (6*20)  must  be 
dimensioned  in  the  calling  routine;-  some  of  the 
elements  of  this  array  are  not  used  In  SMEAR  3] 

Original  density  (g/cc) 

Bulk  modulus*  C  (dyn/cm3) 


D  and  S  In  the  pressure  equation: 

f  -  Cu  ♦  Dw2  ♦  Su3,  where  y  -  DH/ESC(M,1)  -  1 
2 

Shear  modulus  (dyn/cm  ) 

IB 

Gruneisen's  ratio 

Array  contains  number  of  shear  bands,  bsnd  lengths, 
and  plastic  strains  for  each  cell.  The  number  of 
array  elements  required  for  each  cell  is 
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TA&LK  A-l 


(Conclude) 

Vorloblo _ H— nlna _ 

"at  "  3I'i  ♦  »<*i> 

vhoro  Ft  •  0  If  FMUC1  -  0 
»  1  if  FMUCi  >  0 

1  -  0  for  m(N,31)  -  0 

-  I?1  for  m(M,31)  -  1 

-  4  for  »n(M,31)  -  2 
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systam  of  UMP,  the  1.1*4  ot  about  lino  132  should  bo  changed  to  1.1-8. 
Also  tho  Internal  enarg 7  wot  bo  transformed  frooi  erg/g  to  Hbor. 

Po*T80TT  "  *HENP  x  1q1 

Lika  HBff,  EPIC  uooo  tho  energy  por  unit  volume  definition  for  E.  Tho 
input  quantities  ora  do flood  in  Tab la  A-2,  with  nominal  values  in  tho 
units  usod  in  TKSTT 

fur  triggering  a  slide  line  based  on  tho  shear  band  damage  level, 
the  parameters  BF1(M,32)  and  TAD  ore  usod.  TAD  is  intended  to  vary  from 
sore  to  approximately  one,  with  values  of  0.3  or  larger  being  indicative 
of  sufficient  damage  to  causa  full  separation.  The  value  of  TAD  can  be 
associated  with  a  particular  direction  of  damage  (shear  band  orientation) 
or  with  all  damage  i  this  cholca  is  made  through  selection  of  EPR(M,32) . 
For  BFl(M,32)  set  to  I,  where  I  equals  1,  2,  3,  or  4,  the  Xth  shear  baud 
orientation  la  used.  For  a  problem  in  which  a  mlsalle  were  penetrating 
in  the  X-dlrectloo,  gVt(M,32)  could  be  sat  to  1  to  associate  TAV  with 
bands  in  the  plugging  direction. 

The  total  number  of  extra  variables  required  per  cell  will  depend  on 
the  number  of  shear  band  orientations  used  (governed  by  FMUC^),  the  value 
of  BFR(M,31) ,  and  the  variables  already  present  in  the  calling  program. 
With  Ql  in  Table  A-2,  a  formula  la  given  for  the  array  sise  of  CN.  Other 
variables  that  may  not  be  available  normally  are  TAD,  EP,  YHL,  and  ROT. 

F  and  DAOT  are  computed  Just  before  the  call  and  need  not  be  stored. 

2.  Analytical  Form  for  Shear  Band  Site  Distributions 

The  storage  requirements  for  the  shear  band  subroutine,  SHEAR2,  wets 
decreased  by  imposing  an  analytic  form  to  the  sise  distributions.  The 
form  chosen  was  the  exponential  one  that  fits  the  shear  band  data  and 
is  also  like  that  used  for  ductile  and  brittle  fracture.  The  number  of 
bands  per  cubic  cm  with  radii  greater  than  R  is 

-  Nq  exp  (-R/Rj)  (A-l) 
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TABLE  A-2 


XHPUT  DATA  FOB  SHEAR  3 


Nominal 

Variable  _ Values  .  Definition 


% 

BFR(M,22) 

10  (nondla.) 

Cq:  (roach  coefficient  in 
dR/dt  -  CcRdc/dt 

BFB(N,23) 

3 

5  par  ca 

Ajj,  nuclonclon  coefficient 

BFR'M,24) 

0.002  ea 

R  ,  nuclonclon  nix*  paraaoCor 
o 

BFR(M,23) 

2  (nondin.) 

n,  axponanc  In  nuclonclon  race 
•qua  don 

BFR(M,2S) 

0.2  (nondla.) 

C  ,  nuclonclon  threshold  strain 

n 

BF1(M,27) 

0.07  (nondla.) 

b,  ahoar  distortion  assoclstod  with 
a  band 

BFR(H,28) 

0.01  ea 

Rg,  aaxlaua  band  also  at  nucloatlon 

LFR(M,29) 

I  10. 0  (nondla.) 

V^,  daaaga  factor  In 

T  -  V__.b.6irN  R? 

FR  o  1 

BFR(M,30) 

0.377  (nondla.! 

Tan  whora  $  is  the  internal 

friction  angle 

BFR(M,31) 

2  (nondla.) 

Indicator  for  storage  of  plastic 

p 

strains  £. 

0  no  storage 

1  store  strains  for  FNUC^  +  0 

2  store  all  strains 

BFR(M,32) 

► 

1  (nondla.) 

Designator  for  TAU.  For 

1  •  BFR(M,32) ,  TAU  -  T^i) 

FNUC 

1.0  (nondin.) 

i 

Indicator  array  and  weighting 
factor  for  orientations  to  be  used 
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where  N  is  ths  total  nua»>er/ cm3  and  R.  la  tha  alas  parameter  of  tha 
o  ^ 

distribution.  Tha  bands  ara  presumed  to  ba  circular,  with  a  radius  R, 
and  to  be  in  a  place.  As  in  \e  aarliar  version  of  SHEAR2,  the  bands 
are  assumed  to  nucleate  and  grow  in  aix  discrete  planes,  thus  representing 
approximately  a  complete  angular  distribution.  This  analytic  size  distri¬ 
bution  in  Eq.  (A-l)  replaces  the  representation  by  a  series  of  size  groups 
in  each  angular  orientation.  Each  orientation  direction  is  now  rapreaented 
by  three  numbers,  No,  ,  and  Rg,  where  Rg  is  a  maximum  size  of  the  bands. 

The  following  paragraphs  describe  the  procedure  used  for  nucleation 
of  a  new  distribution,  growth  of  an  existing  distribution,  and  cood>lnation 
of  the  nucleated  and  growing  distribution  to  a  single  exponential  in  the 
form  of  Eq.  (A-l). 


Shear  band  nucleation  is  described  by 


dt 


(A-2) 


PS 

where  N(,^,  de^g  are  number  of  bands  and  plastic  shear  strain  in  the  $0 
angular  direction,  and  Aq  and  n  are  constants.  A  limit  on  the  rate  of 
nucleation  is  provided  by  requiring  that  the  total  distortion  represented 
by  the  new  bands  not  exceed  the  total  imposed  plastic  snear  strain.  The 
shear  strain  associated  with  nucleation  is 


.  PS  TT  f  bR3dN 

mjL 

where  b  ■  B/R  is  the  maximum  relative  offset  across  the  band 
the  number  nucleated  with  a  radius  R.  By  using  Eq.  (A-l)  to 
derivative  of  N  in  terms  of  R,  and  integrating  Eq.  (A-3),  we 


(A-3) 

and  dN  is 
obtain  the 
obtain 


AePS  -  3ffbANR3 
n 


(A-4) 


where  AN  is  the.  ‘  otal  number  per  cubic  cm  nucleated  with  a  size  distri¬ 
bution  parameter  R^  ■  R^.  Eq.  (A-4),  is  then  an  upper  bound  on  the  number 
nucleated  in  any  time  interval  . 


t 
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Growth  of  the  shear  bands  la  given  by  tha  viscous  law  from  Reference  18. 


dcPS 

C  a  _ii 

l‘<ri  dt 


(A-5) 


where  la  a  growth  coefficient.  When  this  viscous  growth  process  is 
applied  to  a  else  distribution  as  in  Rq.  (A-l),  the  entire  distribution 
Increases  in  else.  The  resulting  distribution  has  the  sane  exponential 
fora,No  is  unaffected,  and  the  change  In  R^  is  obtained  through  the  use 
of  Xq.  (A-5). 

Two  constraints  are  applied  to  the  growth  law  to  Halt  the  maximum 
growth  velocity  to  the  sound  speed  and  to  require  that  the  total  distortion 
represented  by  the  growth  does  not  exceed  the  total  Imposed  plastic  shear 
strain.  The  Maximum  growth  rate  occurs  at  the  maximum  radius,  according 
to  the  viscous  growth  law,  and  that  radius  is  called  R^.  Hence,  the 
velocity  limit  implies  that  the  Maximus  radius  at  the  end  of  the 
Interval,  Rg2, 


*k  i  Si +  V* 

where  Is  the  shear  wave  velocity.  The  increase  in  the  radius  parameter 
Is  limited  by  maintaining  a  constant  ratio  of  R^/R^.  The  distortion 
limit  cn  growth  is  derived  in  the  same  manner  as  in  Equations  (A-3)  and 
(A-4).  The  Increase  In  plastic  strain  associated  with  the  band  growth 
on  the  00  i>lane  is 


“♦9  -  '  RW> 


(A-7) 


where  R2  is  the  final  size  distribution  parameter.  Hence,  based  on  the 
imposed  distortion,  this  limit  is 


r: 


ArPS 


106  3TTbN 


006 


(A-8) 
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iAa 


Casblaiag  Equations  (A-5),  (A-4) ,  sad  (A-0),  m  find  ths  vslus  of  Rj 
on  ths  #6  plsas  to  bo 

*»♦*-  [(*M9  +  V‘>  • 

*i*e  •xp(ccASe)  *  (A“9) 


Aftsr  calculating  growth  of  ths  old  also  distribution  and  nudoation 

of  a  asw  distribution  with  slsa  paraaator  1  »  wo  conblno  tho  two  diatri- 

n 

butlons  to  fona  a  now  distribution  with  tho  sons  form  as  Eq.  (A-l).  Tho 
cosbinsd  distribution  hss  tho  paranotors  N,.  and  Ry  Tho  cosiblnatloa  is 
sndo  with  two  constraints: 

o  Ths  now  curvo  prosonts  tho  also  and  auaber  of  tho 
largost  shoar  bands  (K  «  1^). 

s  Tho  now  curvo  roprosonts  tho  sun  of  tho  plastic  strain 
of  tha  nudoatod  and  growing  distributions. 

Tho  first  roqulroaont  statoe  that  tho  nunbar  of  bands  at  is  prasorved; 
honco 

M3B  "  V^VV 

(A-10) 

»  N^xpC-^/Rj)  +  ANoxp(-RB/Rn) 

Tho  total  plastic  strain  of  tha  now  distribution  is 

ePS  -  3rbN3R^  (A-ll) 

By  eliminating  N3  fron  Eq.  (A-ll)  with  the  aid  of  Eq.  (A-10),  we  obtain 
an  expression  for  Rj. 
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(A-l  2) 


ePS  -  3irbN3Bexp  (R^/R^)  R* 

where  Hjb  Is  known  from  Zq.  (A- 10) .  A  staple  algorithm  for  solving 
Eq.  (A-12)  was  developed  for  the  progrsn.  N,  is  then  obtained  from 
Eq.  (A-ll) ,  and  the  final  size  distribution  is  fully  specified. 

The  foregoing  nucleation,  growth,  and  combination  procedures  have 
been  Incorporated  into  a  version  of  the  shear  band  routine  SHEAR2.  We 
have  made  a  series  of  calculations  and  are  comparing  then  with  results 
froa  the  earlier  version.  The  Initial  observations  show  that  there  are 
differences  in  the  results  of  the  two  versions.  These  differences  stem 
primarily  from  the  approximation  (in  the  new  routine),  involved  in  the 
recombination  calculations  above,  and  the  summation  approximation 

ZNiRi  "  J  r3<IN  (A-13) 

implied  in  the  original  version.  The  sum  taken  over  all  the  size  groups 
is  a  poor  estimate  of  the  Integral  unless  a  large  number  of  groups 
(at  least  10)  are  used  and  are  spaced  along  the  function  in  an  optimum 
manner. 

3.  Stress  Relaxation  Algorithm 

The  new  SHEAR3  stress  relaxation  process  caused  by  a  combination  of 
yielding  and  growth  of  shear  band  damage  is  derived  here.  This  process  is 
part  of  a  complete  computational  description  of  the  stress-strain  relations 
for  a  material  undergoing  shear  band  damage.  The  basic  assumptions  of 
the  model  arc  presented  first,  followed  by  yield  calculations  for  one 
or  more  slip  planes.  The  possibility  of  a  combination  of  shear  and 
tensile  damage  is  treated.  In  each  case  the  relaxed  stresses  and  plastic 
strains  are  obtained.  Then  steps  are  provided  for  implementing  the 
stress  relaxation  algorithm  into  a  shear  band  model. 


Hh  a«J«t  of  tbs  itrui  relaxation  process  is  that  the 

shear  strength  of  a  partially  shear  banded  notarial  depends  jointly  cm 
the  amount  of  shear  banding  and  the  nernal  stress  across  the  bands. 
Quantitatively ,  the  shear  strength  is 

» •  s  °D  \ 

where  i  is  the  band  length 

h  is  the  length  of  the  notarial  block  in  the  band  direction 

a  is  the  nomal  stress  across  the  band 
H 

d  is  the  friction  angle 

Y  is  the  yield  strength  in  shear  (one-half  the  tensile 
yield  strength). 

The  configuration  of  the  band  and  stresses  is  shown  in  Figure  A-l. 

Wots  that  the  band  reduces  the  shear  stresses  S12  and  S21  equally  but 
has  no  affect  on  the  *23  stresses. 

Equation  (Arid)  is  extended  to  the  shear  band  else  distributions 
by  Introducing  the  danage  density  function  T^,  which  represents  the 
total  area  of  shear  bands  per  unit  voluno  associated  with  the  k-th  plane. 
The  notarial  is  fully  shear  banded  for  Tk  ■  1.  Then  Eq.  (A-1A)  is 
assumed  to  take  the  fom 

^  -  xk  0atan  d  ♦  (1  -  Tk)Y#  ,  ^  <  Y#  (A-15) 

The  shear  strength  is  always  positive  and  cannot  exceed  Y  . 

8 

Central  to  the  use  of  Eq.  (A-15)  are  the  assumptions  t  t 

o  The  shear  band  danage  can  be  allocated  to  planes 
of  interest  (such  as  the  X-Y  plane)  Instead  of  to 
specific  directions  (as  in  tensile  fracture). 

o  The  yield  process  can  be  dealt  with  on  specific 
planes  (a  Tresca  yield  approach),  rather  than  by 
considering  all  stresses  simultaneously,  as  in  a 
Mlses  yield  process. 

Both  these  assuaptlons  seen  to  fit  the  basic  mechanisms  mentioned  earlier. 
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Za  tha  yield  procua  it  la  furthar  assuasd  chat  ylaldlng  doas 
aaC  affact  tha  praaaura  vhlla  tha  material  la  undar  compression.  Thu* 
tha  pressure  remains  a  function  only  of  eha  density  and  internal  anargy 
and  la  not  aodlflad  by  ylaldlag  except  by  tha  haat  generated. 

lhaaa  aaayptlooa  ara  tha  basis  of  tha  nodal  darlvad  hara.  With 
thaaa  assumptions,  tha  nodal  will  provide  for  an  anlaotropy  of  strength  so 
that  on  plaaaa  that  ara  shaar  bandad  tha  ahaar  atrangth  nay  go  to  taro, 
but  ahaar  atraaaaa  on  other  pi ana a  will  rataln  full  atrangth.  This  atrass 
raductlon  will  allow  ahaar  bands  to  prepagata  through  conputationai  calls 
repraaantlng  natarlal  undar  high  shaar. 


la  tha  shaar  band  nodal,  danago  in  tha  forn  of  shaar  bands  is 
nuclaatad  and  grown  on  tha  spaclflc  planas  shown  In  Flgura  A~2.  Tha 
ralatlva  shaar  band  araa  on  ona  of  thaaa  planas  la  determined  fro*  tha 
shaar  band  alas  distribution  for  that  oriantatlon 


b.  Determination  of  tha 


•  H0  axp  (-1/lj) 


(A-16 


where  la  the  number  greater  than  R  per  unit  volume 

Nq  la  tha  total  number  par  unit  volume 

B  la  tha  radius  of  a  shaar  band,  assumed  to  be 
circular  and  plana 

R^  Is  a  shape  parameter  of  tha  distribution. 

The  araa  la  than 


A  -  IT 
P 


R2dN 


2ttN  R? 
o  1 


(A-17 


where  dN  Is  the  nuabar  par  unit  volume  with  a  radius  R.  To  nondiraen- 
sionaliaa  the  damage  factor  and  allow  it  to  approach  unity  for  full 
damage,  a  length  factor  c  is  introduced.  Tha  relative  damage  on  the 
plane  is  than 


T  ■  cA 
P  P 

where  T  is  noudimenslonal. 


2ttcN  R? 
o  1 


(A-18 


i] 


Ones  we  ten  tb*  4«tfn  (acton  T  (or  each  (tear  band  orientation 

P 

piano,  a  composite  d— go  factor  la  required  on  the  plena*  of  interaat  for 
the  atrees  calculation*.  The  contribution  fron  each  ahear  band  orientation 
group  to  the  eonpoalt*  factor  la  a  function  of  the  nine  of  the  angle  9 
batuaon  the  aoraal  to  the  plana  of  Interaat  (i-J)  and  the  normal  to  the 
p-th  ahear  band  orientation  plana,  lance,  for  the  i-J  plane 


-  £ 

P-1 


■In  ♦, 


PU 


(A-19) 


The  aun  run*  over  all  four  orientation  group*  uaad  In  the  present  ahear 
hand  nodal  IHIAKl.  With  this  definition  of  the  eonpoalt*  danage, 
for  the  X-Y  plana  receive*  the  rp  value*  for  the  X,  Y,  and  X-Y  directions 
with  a  weighting  factor  of  1.0  and  for  the  XZ  and  YZ  directions  with  a 
factor  of  O.S.  The  snaar  hands  normal  to  the  Z  direction  do  not  contribute. 

The  pianos  of  internet  for  evaluating  the  damage  factors 

are  the  plena*  of  principal  shear  stress.  For  two-dlnanalonal  problems, 

on*  of  these  plana*  la  always  the  X-Y  plana  because  oa  is  always  a 

principal  stress.  (Z  is  taken  to  be  in  the  thickness  direction  far 

planar  problem  and  in  the  circumferential  direction  for  axlsymetric 

problem.)  however,  the  other  two  planes  can  be  at  arbitrary  angles  6 

and  90  +  0  for  the  X  direction  as  shown  in  Figure  A-2.  The  weighting 
2 

factor*  sin  *  «.  are  derived  fron  the  rule*  for  spherical  triangles. 

PXJ 

The  required  weighting  factors  are  given  in  Table  4-3.  With  these  factors 
the  strength  on  any  principal  plane  can  be  evaluated  with  Eq.  (A- 15) . 

c.  Yleldins  on  a  Sinaia  Plane 

The  yield  process  associated  with  Eq.  (A-15)  nay  be  visualized 
with  the  aid  of  a  Mohr  diagran  as  in  Figure  A-3.  The  ordinate  is  the  shear 
stress  or  a  plana,  and  the  abscissa  is  tho  normal  stress.  A  circle 
represents  all  states  of  stress  in  any  orientation  on  a  plane.  The 
center  of  the  cirda  is  at 


o  +  a 
_* _ I 


(A-20) 
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TAILS  A-3 


HEIORIMG  PACTOtS  tin2  *  FOX  KILATIlfG  DAMAGE 
TO  PftXNCIPAL  PUKES 


lild 

Orientation 

tin2  » _ 

2-Z  Plane 

1  X 

2* 
cos  6 

aln26 

1  Y 

1 

2 

•me 

2 

cos  6 

e  i 

0 

l 

1 

2  ZY 

0 

1/2(1  ♦  sin  26) 

1/2(1  -  sin  26) 

4  YZ 

1/2 

1-1/2  coe2e 

1-1/2  ain20 

3  XZ 

1/2 

1-1/2  sin2 8 

1-1/2  coa2e 

Moca:  6  is  th«  angle  to  principal  direction  1  in  the  X-Y  plane, 
neaaured  counter-clockwise  fron  X  toward  the  T  direction. 

*Ho  band*  are  currently  provided  in  the  Z  direction  because 
there  la  no  shear  on  the  plane. 
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MAWM 

FIGURE  A-3  MOHR  DIAGRAM  POM  STRESSES.  INSETS  SHOW  ORIENTATION  OF 
STRESS  AT  A  POINT 
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where  0^  and  in  principal  atraaaaa  that  determine  tha  Intercepts  of 

tha  . Ire la  with  tha  axis,  and  o  and  o  art  any  two  normal  atraaaaa  on 

*  * 

rrtho ponal  planaa.  Tha  email  Insat  blocka  in  Figure  A- 3  ahow  tha 

c~laatatloa  of  atraaaaa  In  tha  positive  algn  convent ion  correapondlng 

to  points  on  tha  Mohr  circle.  A  couaterelockvlaa  rotation  of  6  to  tha 

P 

direction  of  principal  stress  la  represented  by  a  clockwise  rotation  of 
26^  co  the  Mohr  diagram. 

The  yield  limit  Is  represented  by  a  sloping  Una  In  the  Mohr 
diagram  up  to  a  shear  stress  of  Y#  and  a  horlsontal  line  for  larger 
normal  stresa.  The  slope  Is  similar  to  the  usual  Coulomb  friction  law 

8  ■  0„  tanC  ♦  C  (A-21) 

n 

where  5  la  the  angle  of  Internal  friction  and  C  is  the  cohesion. 

Comparing  Iqa.  (A-21)  and  (A- 15) ,  it  Is  seen  that 

tan  C  ■  tan  +  (A-22) 

c  -  (1  -  *tin9  (A-23) 

In  tha  sketch  In  Figure  A-4  tha  Mohr  circle  Is  tangent  to  the  yield 
limit  line,  so  the  material  is  yielding.  Held  Is  actually  occurring 
on  a  plane  with  a  stress  state  defined  at  point  A,  at  an  angle  of 
AS*  +  5/2  from  the  maximum  principal  stress. 

From  the  foregoing  yield  process,  the  stress  states  at  yielding 
can  be  fully  determined.  As  shown  in  Tersaghi,^  the  principal  stresses 
on  the  Mohr  circle  at  yield  are 

-  2C^S^  +  02  (A-24) 

where  -  tan2 (45*  +  5/2)  -  (1  ♦  sin5)/(l  -  sln5)  and  is  called  the 
flow  value.  In  general  then,  for  the  three  p^islble  planes  between  the 
principal  directions. 
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FIGURE  A~4  MOHR  CIRCLE  FOR  STRESS  STATES  WITH  A  YIELD  LINE  FOR 
YIELDING  ON  THE  1-2  FLANE 
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°i  * J(l  -  VT»  V^eZI +  °j  "( 


(A-25) 


where  and  axa  tha  daaaga  factor  and  flow  value  of  tha  1J  plana 
and  >  Oj. 

For  computational  purposes,  It  la  convanlant  to  dafina  a  yield 
function  g^  that  la  saro  at  yield,  negative  for  elaatlc  behavior,  and 
positive  for  inadniaaible  state a  beyond  the  yield  liait: 


*ij  "  “*  [  “  (l  “  TiJ)Ya  VNC ij  “  °J 

0i-01  y 

2  “  » 


<N  -  1) 


(A-26) 


This  definition  of  g  lncludea  the  liait  on  the  yield  strength.  Tha 
value  of  g  represents  tha  excess  shear  strength  above  the  yield  line. 

Having  defined  vhe  yield  function  g^ ,  we  can  determine  the  post¬ 


yield  stresses  and  plastic  strains.  Let  the  principal  ati 


computed 


by  assuming  all  the  strains  are  elastic  be  O  and  the  yield  function 

N  * 

based  on  these  stresses  be  g^  .  The  postyield  stresses  are  o^,  where 


°i  "  °i  +  Aoi 


(A-27  ) 


On  the  ij  plane,  yield  occurs  without  disturbing  the  third  principal 
stress  or  the  pressure.  Therefore- 


o+o  **  oN  +  0N 
^  i  1  J 


Aoa  ■  -  AOj 


(A-28) 


_N  _N 


Replacing  o^  and  Oj  in  Eq.  (A-26)  with  o^  o^,  and  Aoi  from  Eqs.  (A-27) 
and  (A-28),  we  obtain 


AOi  "  ‘  Aa j  "  N?lj2+  1  8ij 


(A-29) 
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This  solution  is  slso  valid  for  the  horizontal  portion  of  the  yield 
limit  if  is  set  to  1.0  for  that  range.  The  plastic  strains  are 
computed  from  the  change  in  deviator  stress  associated  with  the  yielding 


These  plastic  strain  increments  are  consistent  with  a  normality  flow 
rule  applied  to  deviator  stresses  only,  not  to  total  stresses  as  in 
the  usual  case.  This  special  treatment  of  normality  occurs  because 
of  the  requirement  that  the  pressure  be  unaffected  by  yielding. 

The  foregoing  analysis  gives  a  complete  solution  for  stresses 
and  plastic  strains  if  yielding  occurs  on  only  one  principal  plane  and 
tensile  failure  is  not  occurring. 

d.  Yield  Calculation  for  Multiple  Yield  Lines 

Yielding  may  occur  simultaneously  on  the  three  planes  associated 
wth  the  principal  stresses  in  pairs:  1-2,  1-3,  2-3.  Such  a  possibility 
is  illustrated  by  the  Mohr  circle  and  limit  lines  in  Figure  A-5.  Each 
of  the  three  circles  is  tangent  to  its  own  limit  line.  Note  that 
the  yield  limits  are  not  identical  in  all  three  planes  as  t  ey  would 
be  for  an  isotropic  material.  Thus  we  permit  the  shear  strength  to 
reduce  on  a  plane  of  damage,  but  to  remain  high  on  other  planes. 

Another  sketch  for  viewing  the  multiple  yield  process  is  shown 
in  Figure  A-6.  Because  the  pressure  is  held  constant,  yielded  and 
unyielded  principal  stress  states  can  be  depicted  completely  in  this 
plane  for  all  three  principal  stresses.  This  sketch  is  made  for  the 
simple  case  in  which  -  0  for  all  planes.  The  Mises  yield  criterion 
is  also  shown  for  comparison.  Three  stress  states  A,  C,  and  E  are  shown 
outside  the  yield  limits:  they  represent  a^N,  and  a ^  values.  In 

accordance  with  Eq.  (A-28)  for  the  stress  corrections  for  yielding,  the 
stress  points  are  moved  onto  the  yield  lines  at  B,  D,  and  F.  For  example, 
AB  f«  vertical  because  is  held  constant  during  yield  to  the  1-3  line, 
but  point  B  is  also  outside  the  1-2  yield  line,  so  B  is  not  a  valid 
solution.  Instead,  the  solution  is  at  G,  at  the  point  between  the  1-2 
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TENSION 


o,  COMPRESSION 


FIGURE  A-5 


L  V 

°3  °7 

NORMAL  STRESS 

MA-7893-6 

MOHR  DIAGRAM  ILLUSTRATING  SIMULTANEOUS  YIELD  ON  THREE  PLANES 
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MA-7893-7 


stress  STATES  (t  -  T  -  0) 
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«nd  1-3  yield  lines.  Clearly  the  stress  reduction  along  the  line  AG  can 
be  resolved  into  components  parallel  to  AB  and  EF  so  the  stress  reduction 
can  be  regarded  as  a  combination  of  yielding  on  both  the  1-2  and  1-3  lines. 
A  similar  combined  result  occurs  for  yielding  at  C.  The  line  EF  represents 
yielding  cm  a  single  line.  Yielding  at  I  represents  an  inconveniently 
complex  case.  The  single  yield  computation  on  the  2-3  line  leads  to  a 
solution  that  is  inadmissible  because  it  shows  yielding  on  the  1-2  line 
although  I  is  on  the  elastic  side  of  that  line.  For  point  I,  the  correct 
solution  is  at  point  H. 

When  multiple  yield  lines  are  involved*  the  solutions  for 
yielding  on  a  single  line  (Eq.  A-29)  are  obtained  first*  and  the  resulting 
stresses  are  used  to  evaluate  the  g^  values  on  the  other  yield  lines. 

If  the  single  solutions  do  not  represent  valid  stress  states  (all  g^<  0), 
then  a  solution  is  made  for  stress  states  at  a  joint  between  yield  lines. 
The  joint  solution  is  obtained  by  eliminating  o'  from  the  three  g^ 
equations  using  the  definition 

o'  +  o'  +  a'  -  0  (A-31) 


Then  the  g^  values  are  set  to  zero  and  solved  simultaneously  in  pairs 
for  o£  and  o'y  The  results  for  the  1-2,  1-3  joint  are 

-  j*."1  -  t13>  Jv 3S£12  +  <l  -  T12 FF  "U31 


+  P(2N^12  /(Npi2  +  Nri^  +  Njri2  Nn^)  (A-32) 


512  r  "513  512  513 ' 


and 


°2  ‘  - 

N£12 


2(1  -  t12)Y, 


F 


(A-33) 


ri2 
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la  each  cam  the  third  principal  atraaa  la  obtained  fron  tha  definition 

of  praaavra 

Oj  •  3P  -  02  -  0X  (A-34  ) 

For  yielding  at  tha  1-2,  2-3  point,  tha  solution  la 

°I  ■  -  V  \n  +  (1  -  Ti2)  J\n  (l  +  Vs'1 

+  -  »£J3  -  1)J  /(I  +  Nr ?3+  *^12  V)'  (*-35) 

Again  o2  and  o3  are  obtained  fron  Bqa.  (A-33)  and  (A-34).  For  yielding 
at  the  1-3,  2-3  joint,  tha  solution  is 


-  |2Y#[(1  -  t13)  (1  +  n^23^  *  T23^  J\23N^13^ 

+  P  (M513  -  »523  -  X)|  /(I  ♦  h-13  ♦  H523) 


-  o1(l  + 


513 


-)  +  2 


(1  ■  T13)Y. 
HC13 


+  3P 


(A-36) 

(A-37 


and  is  fron  Eq.  (A-34).  Tha  atraaa  solutions  in  £qs.  (A-32) , (-33), 
(-35)- (-37)  -are  predicated  on  tha  aasuaptlon  that  tha  sloping  portion  of 
each  yield  function  is  being  used.  However,  the  correct  solution  for 
tha  horizontal  portion  can  be  obtained  by  setting  *  1  and  -  0 
in  these  equations. 

Following  solution  for  stresses  at  a  joint  between  two  yield 
lines,  the  g^  value  for  the  third  yield  line  is  evaluated  with  the  new 
stresses.  This  test  avoids  use  of  solutions  at  a  joint  such  as  point  K 
in  Figure  A-6. 

For  a  point  such  as  L  in  Figure  A-6,  all  three  yield  lines  are 
exceeded.  Therefore,  three  individual  yield  calculations  are  made  and 
tested  and  then  three  joint  solutions.  The  joint  solutions  at  G  avid  H 
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are  squally  aatlafactory;  therefore,  an  average  of  tha  strassas  at  thaaa 
two  points  Is  taken  as  tha  solution. 

Tha  plastic  strains  for  tha  joint  solution  ara  takan  from 
Iq.  (A- 30) . 

a.  Tanalla  Failura  . 

Tonsils  failura  nay  occur  and  altar  tha  strassas  in  a  manner 
different  frcn  that  caused  by  tha  shear  processes.  Hare  ve  consider 
only  a  simplified  tensile  failure,  not  a  complete  analysis  of  tha 
opening  of  the  bands  in  tension. 

As  tha  stress  state  becomes  more  tensile,  tha  Mohr  circle  (as 
in  Figure  A-4)  novas  to  the  left.  The  yield  calculations  will  shrink 
the  circle,  causing  and  to  approach  each  other.  Finally,  whan  the 
canter  of  the  circle  novas  to  tha  left  of  the  intercept  (B)  of  tha  yield 
line  on  the  normal  stress  axis,  the  stress  solution  froa  Eq.  (A-29)  will 
give  Oj  >  J^.  Such  a  atate  does  not  actually  swat  tha  yield  requirements. 
Therefore,  two  adjustments  are  made: 

e  The  stresses  are  aet  equal  to  tha  value  at  tha  intercept. 

e  The  pressure  is  permitted  to  change. 

The  stresses  at  the  intercept  point  are 

°i  "  °j  “  _  (1  “  T^)  Ys  cos  5  (A- 38) 

This  result  is  obtained  froa  the  geometry  shown  in  Figure  A-4.  The 
third  principal  stress  is  unchanged. 

For  the  multiple  yield  line  solution,  each  stress  involved  in 
yielding  is  tested  to  make  certain  that  it  is  greater  (more  compressive) 
than  the  Intercept  values.  Each  stress  is  involved  in  two  yield  lines 
and  must  therefore  be  compared  with  two  Intercept  values. 

The  foregoing  tensile  process  is  intended  to  represent  a  combina¬ 
tion  of  shear  and  tensile  behavior.  If  only  shear  were  considered,  the 
sloping  Couloefc  friction  line  would  be  modified  to  become  horizontal  at 
S  ■  C  in  the  tensile  range.  Then  there  would  be  no  intercept  on  the 


119 


normal  itrtu  axis  and  no  aaad  for  a  apodal  treatment.  However,  th« 
tanalla  fallura  nodal  outlined  above  was  felt  to  better  represent  reel 
notarial  behavior. 

f .  Inolenantation  in  the  Shear  Bard  Model 

Incorporation  of  the  stress  relaxation  algorithm  into  the 

shear  band  nodal  requires  sons  coordination  of  orientations  used.  The 

stresses  stored  for  the  shear  band  nodal  are  oriented  in  the  fixed , 

external  X,  Y,  Z  space.  However,  the  damage  orientations  are  fixed  to 

the  notarial  and  tharefore  nay  gradually  rotate  about  the  Z  axis  during 

a  no real  computation.  The  stress  and  d ana go  quantities  are  brought 

together  on  the  orientations  of  the  principal  stress  directions.  The 

principal  stresses  are  at  an  angle  8  (positive  counterclockwise)  from 

P 

the  X  direction  in  the  X-Y  plane.  The  naterial  has  rotated  through  an 

angle  p  (also  positive  counterclockwise)  from  the  original  X-Y  orientation. 

Initially,  both  the  stresses  and  damage  quantities  are  transformed  to 

the  principal  stress  coordinates  through  the  rotations  8  end  6  -  p  , 

P  P 

respectively.  Then  the  yield  and  plastic  strain  calculations  above  are 
undertaken.  Plastic  atralna  are  then  transformed  to  the  material 
coordinate  directions  (through  -  9^  ♦  p),  and  the  stresses  are  rotated 
through  -  0p  to  the  external  X-Y  system. 
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APPENDIX  B 


THRESHOLD  INITIATION  CRITERIA  FOR  ADIABATIC  SHEAR  BANDING 


1.  Constitutive  Instability  Theory 


A  combo  explanation  for  tha  Initiation  of  adiabatic  shear  bands 

is  that  they  are  plastic  instabilities  that  occur  whan  the  thvmal 

softening  overwhelms  the  work  hardening,  thus  causing  the  effective 

hardening  modulus  to  become  negative.  Tha  theory  of  such  instabilities 

37  38 

has  been  recently  summarised  by  sice  and  Clifton.  Another  recant 

39 

literature  survey  for  sltear  banding  phenomena  was  given  by  Rogers. 

This  instability  approach,  modified  by  dynamic  loading  effects,  has 
been  applied  to  our  data,  and  is  described  below. 

If  we  assume  that  the  shear  strength  T  of  an  isotropic  material 
is  a  function  only  of  the  equivalent  plastic  strain  EP  and  the  temperature 
T,  then 

t  -  t  <ep,T) 

As  shown  for  4340  steel  (Rc40)  U-7]  in  Figure  B-l,  this  yield  condition 

defines  a  surface,  and  any  specified  load  path  must  lie  on  the  surface. 

Two  extreme  paths  are  the  Isothermal  hardening,  and  constant  strain 

softening,  paths  shown  in  Figure  B-l  as  solid  and  dotted  lines, 

respectively.  An  actual  load  path  will  usually  lie  between  these 

extremes .  An  important  special  case  ia  the  adiabatic  load  path  shown 

as  the  dashed  curve  in  Figure  B-l.  If  the  plastic  work  is  converted 

to  heat,  and  the  heat  la  not  allowed  to  escape,  then  significant  thermal 

softening  can  occur.  Figure  B-2  shows  the  projection  of  the  estimated 

adiabatic  heating  path  on  the  T-ep  plane  for  three  hardnesses  of  4340 
40 

steel.  It  is  seen  that  the  stress  on  an  adiabatic  load  path  begins 
to  drop  at  strains  in  excess  of  about  20Z. 
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LOADING  PATHS: 


FIGURE  0-1  STATIC  TEMPE RATURE-OtFEWDENT  YIELD  SURFACE  FOR  4340  STEEL,  R,43 


FIGURE  B-2  STATIC  AOiARATtC  YIELD  CURVES  FOR  4340  STEEL 


lbs  MttMMtieil  conditions  (or  constitutive  instability  srs 
illustrated  asst  simply  by  considering  the  one-distensions!  case  shown 
in  Figure  1-3.  In  this  esse  a  pure  shear  wave  propagates  in  the 
direction  of  the  Lagranglan  coordinate  h,  carrying  displacement  u  and 
velocity  v  normal  to  h.  The  equations  describing  the  motion  are 


CP  .  .  I  & 
2  9h 


(B-l) 


2?r 


ii  .  .  i  c* 

at  2  c  ah 


(1-2) 


(B-3) 


(B-4) 


2  o 

idiere  pC  ■  dt/dr  and  where  we  have  assumed  that  the  shear  stress 
is  a  function  only  of  the  plastic  strain  ep.  This  will  be  the  case 
whenever  a  particular  load  path  is  specified  because  the  specification 
provides  an  independent  relation  between  tP  and  T,  and  the  T  -  dependence 
can  be  eliminated. 

Combining  (B-l),  and  (B-2),  and  (1-4)  yields 


(iJlV  -  ^ 

\*J/  *  W) 


(B-5) 


We  next  linearise  this  nonlinear  equation  by  considering  a  perturbation 
on  a  static  state  uQ.  That  is, 

u  •  u ♦  d,  C2  »  C2  (Y  ) 
o  o 


Then  (B-5)  becomes 


—2  •  r  <V  *4 

3t2  2  0  3h2 


(B-6) 
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FIGURE  B-3  ONE-OIMENtiGNAL  SHEAR  WAVE  LOAOING 


Which  has  solutions  of  tha  fora 


fi  •  Ra  exp  lik(h  -  -§)] 


V5 


(B-7) 


Khars  k  is  tha  wavs  nuabar  and  1 


•  FT. 


Thus*  if  C  is  raal  and  positive,  tha  periodic  disturbance  is  stable; 

2 

hut  if  C  la  raal  and  negative,  divergent  growth  occurs  and  the  saallett 
wavelengths  are  the  aost  unstable. 

2 

In  aunaary,  shear  bands  will  nucleate  whenever  C  <  0,  where 


Mat -i  [(&),*(#)(&)] 


(B-8) 


As  shown  In  Figure  B-2,  nucleation  will  occur  under  adiabatic  deformation 
In  4340  steel  at  strains  of  20-30Z,  depending  on  the  hardness.  As  we 
discussed  in  Section  III-B,  this  result  Is  In  good  agreement  with  the 
observed  Initiation  strains  In  tha  cylinder  experiments. 


Relation  (B-8)  thus  appears  to  be  a  good  candidate  for  the  desired 
threshold  Initiation  criterion  under  adiabatic  loading  conditions.  The 
remaining  problem  is  to  define  the  conditions  under  which  the  load  can 
be  considered  to  be  suitably  adiabatic. 


Kinematic  Localisation.  One  way  for  adiabatic  deformation  to 
occur  is  for  the  loading  rate  to  be  high  compared  with  both  the  heat 
flow  rate  and  the  rate  at  which  plastic  waves  can  carry  away  the  strain. 

Let  us  first  examine  the  plastic  wave  effect.  For  the  one-dimensional 
situation  of  Figure  B-3,  the  characteristic  directions  have  velocities 
+  C,  and  the  characteristics  carry  constant  values  of  v  +  /C  dep.  When 
the  disturbance  travels  into  initially  undisturbed  material  (a  simple 
wave),  then  v  and  ep  are  constant  along  the  forward  facing  characteristics, 
and 
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(B-9) 

simple  wave 

(B-10) 
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Combination  of  (B-9)  and  (B-10)  with  (B-3)  and  (B-4)  yields 
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-C 

3ep 

dh 

(B-14) 

Equations  (B-3),  (B-4),  and  (B-9)  through  (B-14)  show  that  Imposed 
high  accelerations  and  corresponding  velocity  gradients  cause  high  strain 
rates  and  strain  gradients;  that  is,  rapid  changes  in  particle  velocity 
will  cause  rapid  changes  in  plastic  strain.  The  effect  of  rapid  mech¬ 
anical  loading  is  most  clearly  seen  from  Eq.  (B-4),  rewritten  as 


To  produce  a  given  strain  gradient,  we  must  have  a  given  value  of 

X  O 

— =•  -57.  If  the  hardening  modulus  is  high,  corresponding  to  a  large  C  , 

2C  dt 

then  a  correspondingly  large  acceleration  3v/3t  is  required  to  cause 
this  "kinematic  localization,"  or  "plastic  wave  trapping." 

If  the  hardening  modulus  is  low,  as  when  the  material  approaches 
the  instability  point,  then  only  a  small  acceleration  is  required  to 
maintain  a  given  dep/dh  or  £p, 

A  necessary  condition  for  wave  trapping  is  that  the  flow  is 
divergent;  that  is,  lines  of  constant  strain  must  diverge.  The  flow 
is  always  divergent  for  a  simple  wave  if  the  material  is  such  that 


because  Chan  large  strains  travel  slower  than  snail  strains.  For  non- 
simple  waves  the  divergence  conditions  any  not  always  be  aet  for  such 
Materials  because  the  initial  particle  velocity  and  strain  fields  nay 
be  convergent. 

In  a  slanted  inpact  experiment,  where  the  target  is  Initially  at 
rest,  the  simple  wave  conditions  apply  for  the  shear  wave,  and  Eq.  (B-12) 
can  be  Integrated  across  the  shear  wave  to  yield 

Av  *  2CACP  (B-12’) 


Because  Aep?  ■  0.2  for  4340  steel,  a  tangential  inpact  velocity  of 
about  0.4  C  will  be  required  to  kinematically  localize  the  critical 
strain  at  the  Impact  surface.  Of  course,  C  is  a  function  of  ep  and 
(12')  should  actually  be  written  as 

v 

Av  -  2  /  C(ep)  dep  (B-12”) 

Nonetheless,  it  is  clear  that  a  high  hardening  modulus  (large  C) 
material  will  require  a  higher  impact  velocity  to  cause  localization, 
whereas  a  low  hardening  modulus  material  will  more  easily  cause  thermal 
failure  at  the  Interface. 

Adiabatic  Heating.  The  second  dynamic  factor  that  determines  if 

deformation  is  adiabatic  is  that  the  material  must  be  heated  significantly 

faster  than  it  cools.  To  determine  the  critical  strain  rate  required 

to  maintain  adiabatic  heating,  we  must  compare  the  rate  at  which  heat 

flows  out  of  a  plastically  deforming  region  with  the  rate  of  heating 

caused  by  the  plastic  work.  To  simplify  the  analysis  as  much  as 

33 

possible,  we  use  Seaman's  definition  of  a  pseudovelocity  of  heat  flow. 


dh  jJk 
dt  h 


(B-15) 
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where  h  Is  now  the  width  of  a  hot  zone,  6  la  a  function  of  the  ratio  of 

temperatures  inalde  and  outaida  tha  hot  zona  (Seaman  takes  g  ■  10  as  a 

representative  value)*  and  k  is  now  the  thermal  diffusivity  < about 
2 

0.15  cm  /a  for  steel).  Note  that  the  heat  velocity  is  si owe  for 
larger  values  of  h. 

The  rata  of  cooling  is  thus  caused  by  the  heat  being  distributed 
over  an  ever-increasing  volume  of  material.  That  is. 


2  £  .  6k 

Q  -h  -„2 

cooling 


(B-16) 


where  Q  is  the  specific  heat  energy.  The  rate  of  heat  increase  is 
thus  the  rate 'of  increase  caused  by  the  plastic  work  minus  the  rate  of 
cooling,  or 


I; P  _  Bit 


Q  • 


(B-17) 


This  equation  can  be  compared  to  the  usual  heat  flow  equation: 


Q 


—  ep  +  k 


(B-18) 


Therefore,  to  allow  approximate  adiabatic  conditions,  we  require  that 
the  first  term  on  the  r.h.s.  of  (B-17)  or  (B-18)  greatly  exceed  the 
second  in  magnitude.  Somewhat  arbitrarily,  we  specify  that: 


;P  >  lOSkQ 

<T/(,)h2 


(B-19) 


That  is,  we  require  that  the  heating  rate  be  ten  times  the  cooling  rate. 

Equation  (B-16)  shows  that  the  relative  rate  of  cooling  is  higher 

for  smaller  widths  h  of  the  hot  zone,  corresponding  to  a  larger  negative 
2  2 

value  of  3  Q/dh  in  the  heat  flow  equation  (B-18). 
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To  specify  the  required  strain  rate,  we  «ust  therefore  specify 
2 

the  value  of  Q/h  that  we  wish  to  maintain.  For  4340  steel,  appreciable 

softening  has  occurred  at  T/T  ■  0.5,  or  for  a  heat  energy  of 

m 

Q  -  3.5  x  10*  erg/g  . 

For  the  contained-fragmenting-cy Under  experiments,  significant 
'strain  gradients  occur  over  distances  of  about  0.1-0. 5  cm.  Relation 
(B-19)  then  gives 

ep  >  102  s"1  -  lO^sec”1  (B-20 


where 

.0  -  10 

2 

k  -  0.15  cm  /sec 

9 

Q  •  3.5  x  10  ergs/g 

T  ■  14  kbar 

3 

p  -  7.85  g/cm 

For  oblique  impact  experiments,  where  an  extremely  thin  layer  near 
the  impact  interface  is  loaded  by  the  shear  wave  via  plastic  vave 
trapping,  widths  on  the  order  of  the  grain  size  of  the  materiai.  may  be 
differentially  heated.  If  the  width  is  assumed  to  be  100  pm,  then 
relation  (B-15)  gives 

ep  >  %  105sec_1  U-20’) 

Thus,  Relation  (B-20’)  can  be  combined  with  (B-4’>  to  determine  the 
Impact  acceleration  required  to  attain  adiabatic  loading  condition  in 
an  oblique  Impact  experiment. 

2.  Development  of  a  Simplified  Formula  to  Calculate  the  Critical 
Strain  Required  for  Shear  Band  Initiation 

In  numerical  computations  the  code  could  simply  check  continuously 
if  inequality  (B-8)  is  satisfied  to  determine  if  t-.he  threshold  condition 
for  adiabatic  shear  banding  has  been  met.  In  practice,  it  is  more 
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enlightening  to  um  a  critical  strain  valua  that  hat  been  determined 
baforahand,  becauaa  wa  ara  than  abla  to  saa  claarly  tha  rola  of  the 
Material  propertlea  of  work  hardening  and  thermal  softening  in 
determining  tha  material  resistance  to  shear  banding.  In  tha  following 
paragraphs  wa  derive  simple  formulae  for  estimating  this  critical 
strain,  and  compare  tha  predictions  of  tha  formulae  to  experimental 
data. 


Analysis.  We  first  assume  that  the  yield  strength  (the  effective 
stress  at  yield)  for  most  ductile  materials  can  be  expressed  as: 


t  -  H(ep)F  ^  |  (B-21) 

— n  + 

where  er  is  the  equivalent  plastic  strain,  H  is  a  strain  hardening 


function,  F  is  a  thermal  softening  function,  T  is  the  absolute  temperature, 
and  T  is  the  absolute  melting  temperature.  For  materials  that  harden 
primarily  by  dlalocatlon-dlalocatlon  interactions,  H  can  be  described 
by  a  power  law  of  the  form: 


H(eP)  -  T  £*"  (B-22) 

o 

where  n  is  the  atraln  hardening  exponent  and  To  is  a  constant. 

We  assume  that,  under  quasl-adlabatic  dynamic  loading  conditions, 
all  the  plastic  work  is  converted  into  heat.  Thus,  the  thermal  softening 
function  F  can  be  expressed  as: 

r(i)  ■  T(k)  ;b-23> 


The  effective  stress  is  defined  as  V-!  ls  the  ij 

^component  of  deviator  stress  tensor. 

The  equivalent  plastic  strain  is  defined  as  Jyj\  dep  dep  ,  where  ep 
is  the  ij  component  of  the  plastic  strain  increment.  1 
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where  W  Is  the  specific  energy  of  plastic  deformation  an d  E^  la  the 
specific  Internal  energy  at  Incipient  aelt.  Experlnental  data  has 
shown  that  the  thermal  softening  behavior  of  aeny  aaterlals  can  be 
fitted  either  by  a  parabolic  expression: 


■w-t-tr 

or  by  a  linear  expression: 


(B-24) 


(B-25) 


where  a  and  a"  are  constants.  The  approach  is  the  saae  whether  a 
parabolic  or  linear  thermal  softening  function  is  chosen;  hence,  only 
the  formula  using  a  parabolic  thermal  softening  function  will  be 
developed  here.  Combining  Eqs.  (B-22)  and  (B-24)  gives  an  expression 
for  the  yield  strength: 


The  next 
depends  on  ep 
The  lncreaent 
given  by: 


t  •  t  e*"1 


(■-?) 


1/2 


step  Is  to  determine  how  the  thermal  softening 
so  that  T  can  be  expressed  as  a  function  of  ep 
In  the  specific  energy  of  plastic  deformation. 


(B»26) 

function 
alone . 
dU,  is 


dV 


1°  rpn 
P 


(B-27) 


where  p  is  the  density  of  the  material.  Rewriting  Eq.  (B-27)  gives: 


(B-28) 


Integration  of  the  left-hand  side  of  Equation  (B-28)  between  the  limits 
of  0  and  V,  and  integration  of  the  right-hand  side  between  the  limits 
of  0  and  y  gives: 
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(B-29) 


1  - 


UTt P 
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Tt¥l 


2p  E  (n+1) 


Combining  Eqs.  (B-29)  and  (B-26)  yield*  an  expansion  for  T  in 
tans  of  cp  alone: 


t  •  t  epn 
o 


aT  2-n(2iH.l) 
o 

2pE  (n+1) 


(B-30) 


Equecion  (B-30)  describee  the  effective  stress  at  a  given  point 
in  the  aaterial  in  tens  of  the  equivalent  plastic  strain  at  that 
point  under  quasl-adlabatic  dynamic  loading  conditions.  Equation  (B-30) 
also  describes  the  macroscopic  dynamic  stress-strain  behavior  of  the 
material  as  long  as  the  deformation  remains  homogeneous.  However,  once 
the  macroscopic  thermal  softening  behavior  overwhelms  the  macroscopic 
strain  hardening  behavior  (that  is,  when  <  0),  a  constitutive 
instability  la  triggered  and  the  deformation  becomes  localized. 

In  other  words,  we  assume  that  ahuar  banding  initiates  when  the 
conditions  for  constitutive  instability  are  reached.  As  discussed  in 
Section  III-B,  the  number  of  points  in  the  materia1  at  which  shear  bands 
actually  do  nucleate  is  governed  by  a  separate  nucleatlon  function,  and 
the  relative  growth  of  the  nucleated  shear  bands  is  governed  by  a 
separate  growth  function.  .  In  any  event,  shear  banding  initiates  when: 


dT 


diP 


nT 

o 


o(2n+l)T^ep 

(rH-1)  2pE 

ni 


(B-31) 


The  critical  strain  required  for  the  initiation  of  shear  banding  under 
quasi-abiabatlc  dynamic  loading  conditions  is  therefore  given  by: 


(B-32) 
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If  a  linear  thermal  softening  function,  Equation  (B-25),  la  used, 
cjr  ia  glvan  by: 


Tha  coapatlng  affacta  of  thermal  aoftaning  coefficient,  yield  strength, 

and  strain  hardening  coefficient  are  evident:  Increasing  a  increases 

cp  ,  whereas  Increasing  t  and  a  decrease  ep  . 
cr  o  cr 

Application.  We  have  used  Bq.  (B-32)  to  predict  the  critical 

strain  required  for  shear  band  Initiation  In  five  different  steels. 

Tbs  results  are  surprised  in  Table  B-l .  The  specific  Internal  energy 

at  Incipient  salt  was  taken  to  be  the  enthalpy  of  the  material  at  the 

temperature  In  which  the  material's  strength  has  gone  essentially  to 

aero.  This  turned  out  to  be  about  the  same  for  each  material.  The 

sources  from  which  the  various  values  of  T  ,  n,  and  a  were  extracted 

o 

for  each  material  ara  indicated  at  the  bottom  of  the  table.  The 
materials  in  the  table  are  listed  according  to  their  predicted  resistance 
to  shear  band  initiation  during  ballistic  impact  by  blunt-nosed 
projectiles,  with  R^O  4340  steel  being  the  least  resistant:.  Hots 
from  the  last  column  in  the  table  that  Equation  (B-32)  predicts  that 
Rc40  4340  steel  will  initiate  shear  banding  at  a  relatively  low  value 
of  equivalev  plastic  strain  (M.9Z) .  On  the  other  hand.  Equation  (B-32) 
predicts  that  mild  ateel  will  initiate  shear  banding  only  at  a  very 
high  value  of  equivalent  plastic  strain  (^1381) . 

The  ranking  of  the  relative  resistance  to  shear  banding  predicted 

by  Equation  (B-32)  for  the  materials  listed  in  Table  B-l  agrees  with 

31 

experimental  observations.  Erlich  et.  al.  observed  that;  R  40  4340 

c 

steel  shear  bands  *•'  tantially  at  equivalent  plastic  strains  less 

than  402.  Sriicb  a-  also  observed  that  R  22  RHA  does  not  shear-band 

c 

until  equivalent  plastic  strains  greater  than  about  80Z  have  been  reached 

(except  along  planes  of  weakness  parallel  to  the  rolling  direction)  . 

41 

Finally,  Shockey  et.  observed  that  low-carbon  steel  does  not 

shear-band  even  aft.  extensive  plastic  deformation  under  ballistic 
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Inpact  conditions.  It  thus  appears  that  Equation  (B-32)  can  ba  uaad 
to  correctly  rank  nsterlal  aa  to  their  relative  resistance  to  shear 
band  formation  and  to  roughly  eatlaate  the  critical  strain  required 
for  shear  band  initiation.  It  should  be  noted  that  these  experimental 
"critical"  strains  represent  only  the  strain  at  which  large  shear  banda 
were  first  observed,  and  not  the  critical  strain  at  which  shear  banding 
really  initiated.  If  the  actual  critical  strains  could  be  pinpointed 
•ore  closely  experimentally,  the  critical  strains  predicted  by  Equation 
(B-32)  night  better  natch  experimental  data. 

Conclusions.  A  fornula.  Equation  (B-32),  has  been  developed  for 
quickly  ranking  naterlals  as  to  their  resistance  to  shear  banding  during 
ballistic  impact  and  for  estimating  the  critical  strain  required  for 
shear  band  initiation.  Development  of  this  fornula  required  two  key 
assumptions.  First,  during  quaa 1-ad labs tic  dynamic  loading,  all  the 
energy  of  plastic  deformation  la  converted  into  heat.  Second,  shear 
banding  initiates  in  a  aatarial  when  its  thermal  softening  respore  to 
gnasi-adiabatic  dynamic  loading  overwhelms  its  work-hardening  response 
to  the  same  loading  conditions.  This  fornula  was  used  to  rank  five 
different  steels  as  to  their  relative  resistance  to  shear  banding  under 
ballistic  Impact  by  a  blunt-nosed  projectile.  Thia  relative  ranking 
agreed  with  experimental  observations.  The  estimated  values  of  the 
critical  strain  required  for  initiating  shear  banding  in  each  material 
also  roughly  agreed  with  available  experimental  data..  It  thus  appears 
that  this  formula  has  promise  for  armor  design. 
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4.  How  specifically,  Is  tha  raport  being  usad?  (Information  source,  design 
data,  procedure,  source  of  Idass,  ate.) _ 


S.  Has  tha  in format ion  In  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization,  technical  content,  format,  etc.) 


CURRENT 

ADDRESS 


Name 

Organisation 

Address 


City,  State,  Zip  ’ 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


^ - 

0LD  Organization 

ADDRESS 

Address 


City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


